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ABSTRACT 


The  present  report  describes  the  results  obtained  under  the 
Reentry  Physics  pc’tiun  of  the  Penetration  Systems  Studies  contract. 
The  purpose  of  the  contract  is  to  carry  out  analyses  and  experiments 
of  hypersonic  flow  in  support  of  missile  detection  and  discrimination 
technology.  The  work  comprises  the  following  five  tasks:  (a)  high 
altitude  drag  investigation,  (b)  nonequilibrium  turbulent  boundary  layer, 
(c)  wake  turbulence  measurements,  (d)  improved  laminar  wake,  and 
(e)  near  wake  analysis.  The  last  four  tasks  fall  under  th*»  general 
heading  of  hypersonic  wake  and  boundary  layer  technology.  The  ob¬ 
jective  and  approach  for  each  task  is  presented  together  with  a  summary 
of  the  results  obtained  to  date. 
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SECTION  I 


INTRODUCTION 

Th<“  Reentry  Physics  portion  of  the  Penetration  Systems  Studies 
contract  has  been  a  continuing  analysis  effort  in  the  area  of  hypersonic 
flow  phenomena  with  the  goal  of  supporting  missile  detection  and  dis¬ 
crimination  technology.  The  present  contract  has  five  tasks  outlined 
below. 

One  task  has  centered  around  the  development  of  the  Bird  Monte 
Carlo  technique  to  describe  the  free  molecule/ continuum  transitional 
flow  field  around  hypersonic  bodies  and  carry  out  associated  wind 
tunnel  measurements  to  compare  the  theory  and  experimental  results. 
The  purpose  of  the  task  is  to  be  able  to  obtain  accurate  drag  predictions, 
as  well  as  estimates  of  the  h\gh  altitude  plasma  sheath  properties  around 
a  reentering  body. 

The  remaining  four  tasks  fall  under  the  general  heading  of  hyper¬ 
sonic  wake  and  boundary  layer  technology.  The  second  task  has  been 
the  use  of  a  ballistic  range  to  make  detailed  hot  wire  measurements  of 
tur  bulent  correlation  length  and  spectra  in  the  wake  of  hypersonic 
pellets.  In  this  task,  various  shaped  pellets  are  fired  past  stationa  :y 
single  and  double  hot  wire  probes  to  obtain  the  above  pertinent  statistical 
as  well  as  mean  properties  back  along  the  length  of  the  wake. 

The  third  task  is  to  examine  the  nonequilibrium  turbulent  boundary 
layer  from  the  viewpoint  of  formulating  a  new  low  speed  theoretical 
description,  developing  the  appropriate  low/high  speed  transformations 
and,  finally,  examining  the  effect  of  turbulence  on  chemical  reactions. 
The  fourth  task  has  been  the  prediction  of  a  particular  downstream 
laminar  wake  case  in  order  to  study  the  "breakthrough"  phenomena 
noted  by  L.  N.  Wilson,  General  Motors,^  ^in  a  ballistic  range. 

Finally,  an  analysis  of  the  near  wake  is  being  developed  based 
on  the  Crocco-Lees  critical  point  concept.  This  analysis  considers  the 
upstream  body  boundary  layer  effect  as  well  as  the  pertinent  base  region 
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details.  The  near  wake  work  is  an  outgrowth  of  the  far  wake  turbulence 
analysis  carried  out  during  the  previous  contract  periods.  An  analysis 
of  the  kind  being  developed  is  a  crucial  link  in  the  prediction  of  the  total 
flow  field,  including  downstream  wake,  around  a  hypersonic  vehicle. 

Tht  format  followed  in  the  report  is  to  present  brief  summaries  of 
the  work  where  detailed  reports  have  already  been  issued  during  the 
contract  period.  Where  reports  have  not  yet  been  issued,  or  are  in 
preparation,  as  in  the  case  of  the  High  Altitude  Drag  Investigation,  con¬ 
siderably  more  detail  is  presented  here.  A  list  of  the  previous  reports 
issued  under  the  present  contract  is  given  in  Appendix  I. 


SECTION  II 


HIGH  ALTITUDE  DRAG  INVESTIGATION 

2.  1  INTRODUCTION  AND  BACKGROUND 

Since  radars  now  under  development  are  predicted  to  have  very 
accurate  measurement  capabilities,  the  drag  and  flow  fields  of  reentry 
bodies  at  high  altitude  will  become  important  considerations  in  the  design 
of  future  penetration  aid  systems.  For  example,  the  drag  and  the  nature 
of  the  drag  change  with  altitude  provide  information  about  the  vehicle.  In 
addition,  the  electrons  produced  in  the  stagnation  region  persist  at  high 
altitude  throughout  the  boundary  layer  on  the  body.  This  ionized  layer 
shields  the  aft  corner  thereby  affecting  the  radar  return  markedly  and 
hence,  plays  an  important  part  in  the  discrimination  process.  Thus,  it 
is  important  to  develop  analytical  methods  for  treating  this  transition 
regime  and  to  obtain  additional  experimental  information  about  it. 

The  transitional  regime  between  continuum  and  free  molecule  flow 
has  been  and  continues  to  be  one  that  is  difficult  to  treat  either  experi¬ 
mentally  or  theoretically.  It  is  not  surprising,  then,  that  there  are  so  few 
flow  situations  for  which  a  comparison  between  theory  and  experiment  has 
been  possible.  Most  analyses  have  been  confined  to  flows  with  simple 
boundary  conditions,  for  instance  Couette  and  Rayleigh  flow,  while  the 
few  experiments  deal  with  quite  different  situations  such  as  flows  about 
bodies  and  free  jet  flow.  The  plane  shockwave  is,  of  course,  a  notable 
exception  to  this  situation. 

However,  external  supersonic  flows  about  bodies  are  of  greatest 
practical  interest.  The  objective  of  the  research  reported  here  was  the 
development  of  a  new  theoretical  method  of  treating  such  flows,  the  direct 
Monte  Carlo  simulation  technique  (References  1,  2,  3),  and  a  comparison 
of  theoretical  flow  fields  about  several  simple  shapes  with  those  obtained 
in  a  low  density  wind  tunnel. 

In  addition,  a  comparison  of  the  Monte  Carlo  drag  calculations  on 
cylinders  and  spheres  with  measured  drag  has  been  made. 
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2.  2  DESCRIPTION  OF  THE  MONTE  CARLO  METHOD 

Low  density  supersonic  flow  about  cylinders,  spheres,  sharp 
wedges,  and  sharp  cones  has  been  studied  in  the  free  stream  Knudsen 
number,  Kn  =_jy  range  from  0  025  to  100  and  at  Mach  numbers  from 

3.  8  to  11  using  a  Monte  Carlo  direct  simulation  technique.  This  method 
was  originated  by  Prof.  G.  A.  Bird  of  the  University  of  Sydney.  (Refer¬ 
ences  1,  2,  3)  Prof.  Bird  was  consultant  to  this  project  during  the  past 
year. 


In  the  direct  simulation  technique  a  small  number  of  molecules  is 
set  in  uniform  motion  in  Maxwellian  equilibrium.  A  body  is  instantane¬ 
ously  inserted  in  the  stream  and  the  motion  of  the  molecules  is  computed 
in  a  region  surrounding  the  body.  The  molecular  paths  between  collisions 
are  computed  exactly  but  the  collisions  are  treated  statistically.  The 
calculation  procedure  consists  of  holding  all  molecules  motionless  for  ?, 
time  interval,  At  ,  small  compared  to  a  molecular  mean  collision  time, 
t,  while  collisions  are  computed  everywhere  in  the  field.  Molecules  are 

then  allowed  to  move  with  their  new  velocities  for  the  time  interval  At  , 

m 

and  are  then  held  motionless  in  their  new  positions  while  another  collision 
cycle  takes  place.  Collisions  are  computed  by  statistical  sampling  as 
follows: 

a)  An  estimate  of  the  local  mean  collision  time  t  is  obtained 
from  the  local  density  and  a  characteristic  molecular 
velocity; 

b)  By  comparing  Atm  and  t,  the  appropriate  number  of  col¬ 
lisions  to  occur  during  At  is  known; 

c)  Pairs  of  molecules  are  selected  at  random  from  a  small 
spatial  region  of  extent  -X. (called  a  cell)  and  are  allowed 
to  collide  with  probability  proportional  to  their  relative 
velocity.  This  sampling  scheme  produces  an  average 
distance  between  collisions  equal  to  the  local  mean  free 
path; 

d)  For  pairs  of  molecules  which  have  been  accepted  for  col¬ 
lision,  a  line  of  impact  and  azimuth  in  the  center  of  mass 
frame  of  reference  are  selected  at  random  and  the  con¬ 
servation  laws  are  applied. 

The  initial  flow  field  evolves  to  an  approximate  steady  state  in  this  man¬ 
ner.  Accurate  steady  state  values  of  flow  field  properties  and  body 
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orces  are  determined  by  continuing  this  collision  process  through  a 
arge  number  of  cycles,  computing  a  cumulative  average  of  the  instan- 
neous  values.  In  this  way,  a  small  number  of  molecules  can  be  used 
o  represent  a  large  number  with  accuracy  dependent  only  upon  the  capac¬ 
ity  of  available  computing  machinery  and  the  accuracy  of  toe  molecular 
collision  model  and  surface  reflection  model. 

In  addition  to  the  net  force  on  the  body,  this  method  provides  a  con¬ 
siderably  detailed  description  of  the  flow,  for  example: 

a)  flow  field  properties:  density,  temperature  velocitv 
etc.  ,  everywhere  in  the  region  surrounding  the  bodyf’ 

b)  variation  of  properties  on  the  body  surface-  e  v  ^ 

sure,  skin  friction,  number  flux?  etc.  ;  g*  '  P 

c)  heat  transfer  to  the  body; 

d)  the  distribution  function  of  the  gas. 

A  more  detailed  description  of  the  method  is  given  in  Appendix  II. 
n  essence,  the  method  consist,  of  constraining  a  system  of  a  small  num. 
er  of  molecules  to  behave,  in  the  mean,  as  they  would  if  all  the  other 
molecules  present  in  an  actual  gas  were  present.  'The  Boltzmann  equa- 
mn  is  thus  solved,  in  some  mean  sense,  for  a  rigid  elastic  sphere  gas. 
Effects  of  fluctuations  are  removed  by  continuing  the  calculation  until  a 
arge  sample  of  values  for  any  particular  flow  field  property  has  been 
gathered.  The  mean  of  this  sample  is  then  an  estimate  of  the  flow  field 
property,  to  which  some  uncertainty,  or  variance,  must  be  attached.  A 

more  detatled  examination  of  the  statistical  aspects  of  the  calculation  is 
given  in  Appendix  III. 

In  the  present  work,  an  extensive  development  of  this  method  has 
been  carried  out.  The  theoretical  basis  has  been  further  clarified,  and 
e  computational  techniques  have  been  greatly  improved,  reducing  com¬ 
putational  time  by  an  order  of  magnitude.  Some  innovations  to  the  tech- 
mque  have  been  developed;  for  example: 

?nhe if®1.!  Size\!‘  e*  •  the  size  of  the  small  spatial  regions 

f?eldhl?n  !iaimplJ"g  18  done’  may  now  be  varfed  over  the 
leld,  to  allow  finer  definition  of  regions  of  high  gradi¬ 
ents,  such  as  the  stagnation  region. 
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2)  For  axisymmetric  calculations,  where  the  population 
density  for  molecules  in  our  two-dimensional  compu¬ 
tational  field  must  increase  as  the  distance  from  the 
axis  in  order  to  simulate  axisymmetric  space,  a  sys¬ 
tem  has  been  devised  of  weighting  molecules  propor¬ 
tional  to  their  distance  from  the  axis,  thus  markedly 
reducing  the  required  computer  storage. 

3)  For  two-dimensional  flows,  the  free  molecule  equilib¬ 
rium  temperature  can  be  computed  at  arbitrary  points 
in  the  flow.  Record  is  kept  of  the  energy  flux  and  num¬ 
ber  flux  to  a  small  circle  in  the  flow,  from  which  the 
free  molecule  recovery  temperature  can  be  computed, 
for  an  assumed  type  of  reflection  process.  Since 
molecules  pass  through  the  circle  and  do  not  collide 
with  it,  there  is  no  disturbance  to  the  flow. 

2.  3  MONTE  CARLO  CALCULATIONS 

2.3.1  Tabulation  of  Results 

Digital  computer  programs  have  been  written  to  compute  the  flow 
about  five  simple  shapes: 

flat  plate 

cylinder 

sphere 

sharp  wedge  of  arbitrary  apex  angle 
sharp  cone  of  arbitrary  apex  angle 

A  list  of  the  computer  results  obtained  to  date  is  shown  in  Table  I. 
(No  results  are  presented  for  the  flat  plate. )  All  of  these  results  are 
quantitatively  accurate  to  a  certain  degree.  However,  the  degree  varies 
somewhat  and  for  that  reason  some  are  considered  preliminary,  meaning 
that  either  the  sample  size  is  smaller  than  desirable  or  the  effect  of  some 
parameter  such  as  field  size  or  the  number  of  molecules  in  each  cell 
should  be  investigated  further.  Most  of  the  results  marked  preliminary 
were  executed  by  Professor  Bird  on  the  IBM  7040  computer  in  Sydney. 

The  computer  used  for  most  of  the  other  results  was  the  CDC  6600.  Two 
of  these  results,  the  cylinder  calculation  at  Kn  =  1.  0  and  0.3,  were 
investigated  extensively  to  determine  the  magnitude  of  statistical  scatter 
and  the  effect  of  the  previously  mentioned  parameters  such  as  field  size 
and  movement  time  on  the  results.  This  matter  is  discussed  in  greater 
detail  in  Appendix  III;  however,  it  can  be  said  in  general  that  within  the 
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Table  I.  Computer  Results 


Knud sen 
Number 

Speed 

Ratio 

Body 

Temp. 

Cylinder 

100 

5 

1.  0 

5 

*00 

T+ 

0.  3 

5 

*cu 

Tt 

0.  2 

5 

CO 

Tt 

Preliminary 

0.  1 

10 

• 25  Tt_ 

Preliminary 

0.  025 

5 

tC 

Preliminary 

0.  025 

10 

.49  T 

ten 

Preliminary 

$ 

Tt 

Hi 

5 

xco 

Tt 

igj^ 

■H 

3.  5 

CO 

Tt 

CD 

Preliminary 

1 

5 

■Kin 

un 

5 - 

range  investigated  the  results  are  not  very  sensitive  to  changes  infield 
size,  movement  time  and  number  of  molecules  per  cell.  The  statistical 
scatter  decreases  with  the  square  root  of  the  sample  size;  e.  g.  , 

<rc  *  -pL  .  r 
CD  SSS  D 

where  is  the  standard  deviation.  A  list  of  the  sample  size  and 
standard  deviation  of  the  drag  coefficient  for  the  cases  in  Table  I  is  given 
*n  Table  II.  The  largest  value  of  o^^/C]-}  there  is  2.  6%.  In  general, 
for  a  given  amount  of  computation  time,  drag  coefficient  is  the  most 
accurately  determined  property  of  the  flow  field,  since  it  is  a  gross  fea¬ 
ture  of  the  flow.  Details  of  the  flow  such  as  the  distribution  function  at 
a  point  may  still  show  appreciable  scatter  when  the  drag  is  quite  well 
converged. 

In  the  following  sections  some  of  the  results  in  Table  II  will  be 
examined  in  detail. 

2.  3.  2  Distribution  Function  Calculations 

Liepmann  (Reference  11)  in  a  recent  study  of  the  Bhatnager-Gross- 
Krook  model  applied  to  a  plane  shock  wave  concluded  that  "the  distribu¬ 
tion  function  within  the  layer  is  bi-modal,  exhibits  the  gradual  change 


Table  II.  Sample  Size  Standard  Deviation  of  Drag  Coefficient 


K 

n 

Sample 

r 

-—■(%) 

S 

Tb 

Size 

D 

CD 

Remarks 

Cylinder 

100.  0 

5 

Tt 

2,  193 

2.99 

2.  1 

1 

run 

1.  0 

5 

Tt 

4,  960 

2.44 

1.4 

3 

runs 

0.  3 

5 

Lco 

V 

31,381 

1.85 

0.  5 

4 

runs 

0.  2 

5 

''GO 

Tt 

2,  275 

1.  78 

2.  1 

2 

runs 

0.  1 

10 

.25 

Lao 

Tt 

3,  090 

1. 64 

1.8 

1 

run 

0.  025 

5 

Lco 

Tt 

1, 878 

1.  57 

2.3 

1 

run 

0.  025 

10 

.49 

00 

Tt 

CO 

2,595 

1.  57 

2.  0 

1 

run 

Sphere 

100.  0 

5 

T~ 

.1, 602 

2.  82" 

2.  5 

1 

run 

0.  26 

5 

V 

1,965 

1.88 

2.3 

1 

run 

0.  03 

3.  5 

*'00 

Ttm 

4,769 

1.41 

1.5 

3 

runs 

20u  Sharp  Wedge 

0.  3 

5 

.  09 

Tt 

Lcn 

5,669 

1.67 

1.3 

1 

run 

22.  5J  Sharp  Cone 

0.  1 

5 

L00 

1.47y~ 

1.29 

2.  6 

1 

run 

from  the  molecular  beamlike  behavior  ahead  to  a  Maxwellian  distribution 

behind  the  shock.  The  effect  of  the  fast  molecules  is  noticeable  even 

many  free  paths  behind  the  shock.  "  Liepmann's  results  are  shown  in 

Figure  1,  where  the  distribution  function  f(u,Vv2  +  w  )  is  plotted  versus 

thermal  velocity  u/v  for  several  locations  in  the  shock  wave.  As  we 

ITlflo  # 

proceed  into  the  shock  wave,  f(u)  departs  from  Maxwellian  shape,  the 

peak  shifting  to  the  right  to  lie  near  an  absolute  velocity  equal  to  the  free 

stream  velocity;  e.  g.  ,  the  peak  at  x  =  +2  lies  at  u/vmw=  6.  3,  adding 

U/v  =  3.  0  yields  an  absolute  velocity  U/vm  =  9.  3  which  is  nearly  the 

freJ^stream  Mach  number.  The  sharp  peak  at°positive  thermal  velocity 

thus  represents  a  cluster  of  molecules  which  haven't  suffered  enough 

collisions  to  lose  their  memory  of  the  free  stream.  Deep  in  the  shock 

layer  the  bi-modal  shape  appears.  The  same  qualitative  behavior  of  the 

distribution  function  is  evident  in  the  Monte  Carlo  calculations  of  the  shock 

layer  in  front  of  a  body.  In  Figure  2  is  shown  the  distribution  function 

f(u)  at  several  locations  along  the  stagnation  line  in  front  of  a  sphere  at 

K  =  0.  26,  S  =  5.  In  the  free  stream  f(u)  is  very  close  to  Maxwellian; 

and  as  we  proceed  toward  the  body,  the  bi-modal  shape  appears. 
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Figure  '•  "nr Shock  wav,;  Krook  Modei 

In  Figure  3  comparison  is  made  between  the  calculated  distribute 
function  in  the  uniform  stream  and  a  Maxwellian  at  the  f-,e  stream  tern 
perature,  again  with  very  close  agreement.  The  Monte  Carlo  collision 

calculation  does  produce  a  Maxwellian  distribution  when  the  flow  is  in 
equilibrium  as  it  should. 

2,3,3  Urag  Convergence  with  Sample  Size 

1  Some  examples  showing  the  way  in  which  the  drag  coefficient  con¬ 
verges  are  shown  in  Figures  4  and  5  for  a  cylinder  and  sharp  cone.  Th< 
sample  size  required  for  convergence  to  a  given  level  of  accuracy  will 
vary  with  Knudsen  number  and  body  shape.  Bluff  shapes  and  low  Knudse 
numbers  require  a  larger  sample  than  slender  shapes  and  high  Knudsen 
numbers,  because  bluff  bodies  and  low  Kn  cause  higher  gradients  in  the 
flow,  hence  more  dispersion  from  the  molecular  distribution  in  the  free 
stream.  In  Figure  4  the  results  of  several  calculations  with  different 
values  of  field  size  and  initial  number  of  molecules  per  cell  are  shown. 
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Figure  2.  Monte  Carlo  Calculation 
of  the  Distribution  Func¬ 
tion  on  Stagnation  Line  of 
a  Sphere 
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Figure  3.  Comparison  of  Calculated  Equilibrium  Distribution  Function 
with  a  Maxwellian 

A  similar  but  more  extensive  investigation  for  a  cylinder  at  Kn  =  0.  3  is 
shown  in  Figure  6  and  in  Figure  7  where  one  of  the  cases  of  Figure  6  is 
continued  out  to  an  extremely  large  sample  size.  At  the  end  of  the  com¬ 
puter  run  in  Figure  7  a  total  of  600,  000  molecular  collisions  in  the  field 
had  been  computed.  The  variation  of  the  parameters  in  these  comparisons 
is  not  completely  systematic  and  it  is  impossible  to  completely  isolate 
the  effect  of  any  one  parameter.  However,  it  seems  quite  clear  that 
within  the  range  of  variation  we  have  investigated  their  effect  is  not  large 
and  further  that  after  something  like  2000  to  4000  molecules  have  struck 
the  body  the  drag  estimate  has  converged  enough  for  practical  purposes. 
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Figure  4.  Cylinder  Drag  Coefficient  vs  Sample  Size 
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Figure  5. 


22.  5°  Sharp 


Cone  Drag  Coefficient  vs  Sample  Size 
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Figure  6.  Cylinder  Drag  Coefficient  vs  Sample  Size 
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Figure  7.  Cylinder  Drag  Coefficient  vs  Sample  Size,  Case  17  Continued 


2,3,4  Some  Typical  Flow  Field  Results 

A  collection  of  flow  field  results  is  shown  in  Figure  8  through  13  for 
a  cylinder,  sphere,  sharp  wedge  and  sharp  cone. 

The  cylinder  calculation  is  at  Kn  =  0.  3,S  =  5,  T^  =  T  .  The 
arrangement  of  the  spatial  cells  is  shown  in  Figure  8  and  constant  density 
and  temperature  contours  are  shown  in  Figure  9  and  10.  The  Reynolds 
number  for  this  case  is  27,  based  on  free  stream  conditions  and  cylinder 
diameter.  The  shock  layer  is  a  thick  viscous  region  extending  something 
like  8^  in  front  of  the  body.  Shock  layer  profiles  near  the  stagnation  line 
for  this  and  other  cases  will  be  shown  in  the  next  section.  There  is  a 
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Figure  8.  Cell  Arrangement 
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Figure  9.  Density  Contours  in  Flow  Field  of  a  Cylinder 


Figure  10.  Temperature  Contours  in  the  Flow  Field  of  a  Cylinder 
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large  lateral  extent  to  the  disturbance  ind  a  very  low  density  wake  behind 
the  body  spreading  initially  at  something  like  20  . 

Density  contours  about  a  sphere  at  Kn  =  0.  26,  s  -  5,  Tb  -  Tt^, 
Figure  11,  show  similar  features  with  the  upstream  and  lateral  extent  of 
the  disturbance  reduced,  as  would  be  expected. 

Density  contours  about  a  20°  sharp  wedge  cooled  to  Tb  =  0.  09  Tt^ 
and  a  22.  5°  sharp  cone  held  at  Tb  =  are  shown  in  Figure  12  and  13. 

2.  3.  5  Some  Typical  Shock  Layer  Profiles 

Shock  layer  profiles  of  density,  temperature,  velocity  and  total 
enthalpy  are  presented  for  cylinders  and  spheres  over  Knuds en  number 


range  0.  025  to  0.  3.  The  general  structure  of  the  layer  agrees  with  what 
one  would  intuitively  expect  in  that  a  thick,  viscous  shock  wave  is  present 
and,  for  a  cooled  body,  a  thermal  boundary  layer  exists. 

In  Figures  14  through  18  the  variation  of  density,  temperature, 

velocity  and  total  enthalpy  through  the  shock  layer  on  the  stagnation  line 

is  shown  for  a  cylinder  at  K  =  0.  3  and  0.  025  and  a  sphere  at  K  =0.  26. 

n  r  n 

Temperature  appears  to  be  a  more  sensitive  property  to  detect  the 
beginning  of  the  disturbance.  The  temperature  shock  layer  is*>. 9^  thick 
for  the  cylinder  at  =  0.  3  (Figure  14)  and  ~  6\m  thick  for  the  sphere  at 
Kn  =  0.  26  (Figure  15). 


Two  results  are  presented  at  Kn  =  0.  025  for  a  cylinder,  one  at  s  =  5 
with  T^  =  T^  ,  Figure  17,  and  another  at  s  =  10  with  T^  =  0.  49  T^  ,  Figure 
18.  The  hot  cylinder  temperature  profile  has  a  plateau  at  the  free  stream 


total  temp  T^/Tg,  -  11  which  is  reached  at  about  half  way  through  the 

shock  layer.  The  Reynolds  number  R„  was  364.  This  result  still 

Deo 

shows  an  appreciable  amount  of  scatter  and  is  considered  preliminary. 
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Figure  15.  Sphere  Stagnation  Line  Profiles 


For  the  cooled  cylinder,  the  stagnation  line  temperature  rises  to  the 
free  stream  total  temperature  T^/T^  =11  about  half  way  thru  the  shock 
layer  and  then  falls  in  a  thermal  boundary  layer  to  the  body  temperature. 

=  708  for  this  case.  This  result  shows  considerable  scatter.  For 

eDoo 

both  results  the  velocity  profile  resembles  that  for  a  completely  merged 
viscous  layer. 

For  all  these  results  the  total  enthalpy  risas  through  the  shock  layer 
to  something  like  10%  above  the  f^ee  stream  value  near  the  body.  This 
result  is  still  being  studied. 

The  stagnation  line  properties  for  a  sphere  at  Kn  =  0,03.  S  =  3.  5 
and  T,  =  T.  are  shown  later  in  Figure  34. 
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2.  4  LOW  DENSITY  WIND  TUNNEL  EXPERIMENTS 


2.  4.  1  Sphere  and  Cone  Drag  Measurements 


The  purpose  of  these  experiments,  performed  in  the  TRW  Low 
Density  Wind  Tunnel,  was  to  provide  a  small  set  of  independent  measure¬ 
ments  of  cone  drag  at  Reynolds  numbers  of  thd  order  of  a  hundred.  The 
sphere  drag  measurements  serve,  by  comparison  with  the  extensive  exist¬ 
ing  literature,  to  indicate  that  the  measurement  system  is  accurate  and 
that  the  sting  drag  can  be  properly  accounted  for.  The  wind  tunnel  is 
described  in  Appendix  IV  where  the  method  of  data  reduction  is  also 
given.  The  comparison  of  the  sphere  drag  measurements  with  the  experi¬ 
mental  and  theoretical  data  is  shown  later  in  Figure  36.  The  agreement 
indicates  that  the  present  system  is  satisfactory. 

The  cone  models  had  9°  half  angles  and  were  slightly  blunted.  The 
data  is  shown  in  Figure  19  where  it  is  compared  with  that  of  Potter, 
Reference  19.  The  agreement  at  the  lower  Mach  number,  4.  3,  is  excel¬ 
lent  but  there  is  a  small,  systematic,  and  as  yet  unexplained,  difference 
at  Mach  number  6„  3. 
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Figure  19.  Drag  Coefficients  of  Cones  in  Transition  Flow 
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2.  4.  2  Two-Dimensional  Flow  Field  Surveys 
2.  4.  2.  1  The  Flow  About  Cylinders 


In  the  experiments,  also  carried  out  in  the  Low  Density  Wind 
Tunnel,  flow  fields  in  the  stagnation  regions  of  cylinders  in  argon  and 
nitrogen  flow  were  surveyed  with  the  free  molecule  equilibrium  tempera¬ 
ture  probe.  This  probe,  described  in  References  4  and  5,  is  a  wire,  the 
diameter  of  which  is  small  compared  to  the  mean  free  path,  held  traverse 
to  the  flow.  The  probe  equilibrium  temperature  under  \hese  conditions  is 
a  useful  property  of  the  flow  in  that  it  changes  significantly  through  a 
shock  layer  and  may  be  measured  quite  accurately.  The  probe  is 
described  more  fully  in  Appendix  IV.  The  variation  of  probe  temperature 
with  Mach  number  in  a  uniform  stream  is  shown  in  Figure  20. 

To  provide  information  about  the  character  and  extent  of  disturb¬ 
ances  in  front  of  cylinders  at  various  Knudsen  numbers,  probe  surveys 
were  made  in  the  stagnation  regions  of  cylinders  of  several  sizes.  In 
Figures  21  through  24  lines  of  constant  wire  temperature  are  drawn  to 
scale  with  the  distance  measured  in  free  stream  mean  free  paths.  The 


Figure  20.  Variation  of  Free  Molecule  Equilibrium  Temperature  with 
Mach  Number  in  a  Uniform  Flow  (from  Sherman,  Ref  13) 
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Figure  21.  Lines  of  Constant  Wire  Temperature  in  the  Stagnation  Region 
of  an  Adiabatic  Cylinder  at  =5.4  and  Kn00  =  0.70,  in 
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Figure  22.  Lines  of  Constant  Wire  Temperature  in  the  Stagnation  Region 
of  an  Adiabatic  Cylinder  at  M  =5.4  and  K  =  0.  106,  in 


ifl|lliSii5!«iiS«ii  SliWli  liil  ii 

SmSSffli  *; 

4MHEH8SSSII 


mm  MPS 

fians  iif 


FJiifS  iiHiS!HS  ttlHnHffifl 

liiiil 

:  m  “ 


Figure  24.  Lines  of  Constant  Wire  Temperature  in  the  Stagnation  Region 
of  a  Cylinder  at  =6.3  and  K  m  =  0.012,  in  Nitrogen 
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wire  isotherm  farthest  away  from  the  body  (labelled  A)  is  the  beginning  of 
the  disturbance,  i.  e.  ,  the  location  of  the  first  noticeable  deviation  from 
the  free  stream  values.  Because  the  change  from  the  free  stream  is 
gradual,  the  location  of  A  is  the  least  accurate  of  the  isotherms.  In  Fig¬ 
ure  21  the  body  is  comparable  in  diameter  to  the  mean  free  path  and  the 
disturbance  extends  nine  mean  free  paths  in  front  of  the  body.  Looking  at 
the  values  of  the  temperature  ratio  isotherms,  one  can  see  that  the  shock 
is  just  beginning  to  develop. 

In  Figures  22  and  23,  a  body  at  approximately  the  same  Knudsen 
number  and  Mach  number  is  compared  for  argon  and  nitrogen.  The  extent 
of  the  disturbance  on  the  stagnation  line  is  greater  for  argon  than  for 
nitrogen  while  the  lateral  extent  is  less  for  argon.  It  should  be  mentioned 
that  the  argon  flow  is  strongly  supersaturated;  so  that  while  there  was  no 
indication  of  a  condensation  shock,  it  is  not  possible  to  be  certain  that 
condensation  wa  i  completely  absent. 

Figure  24  shows  the  wire  temperature  contours  about  a  bod/  much 
larger  than  the  mean  free  path;  i.  e.  ,  the  flow  here  is  closer  to  continuum 
flow.  Here  the  extent  of  the  disturbance  and  the  low  values  reached  by 
the  wire  temperature  near  the  body  shows  that  a  shockwave  has  formed. 

It  is  interesting  to  see  that  the  outer  edge  of  the  layer  has  moved  away 
from  the  body,  but  that  the  shape  of  this  disturbance  has  not  changed  much 
xrom  the  example  in  Figure  23.  Hence,  we  start  to  have  a  shock  wave  at 
the  conditions  of  Figure  23  and  at  the  lower  Knudsen  number  of  Figure  24 
we  begin  to  see  a  fully  developed  shock  wave  and  shock  layer. 

Figure  25  shows  the  variation  of  the  stagnation  line  profiles  for 
cylinders  at  various  Knudsen  numbers  in  argon.  The  profiles  display  the 
range  of  shock  layer  structure  from  a  fully  developed  shock  at  Kn  = 

0.  090  to  just  the  beginning  of  a  disturbance  at  Kn  =  0.  70. 

GO 

The  theoretical  value  for  the  recovery  temperature  ratio  of  the 
probe  in  a  uniform  stream  at  high  Mach  number  for  a  monatomic  gas  is 
^w^TT  "  1  =  0.  25  (Figure  20).  In  general,  the  measured  recovery  tem¬ 
perature  does  not  take  on  this  value  in  the  free  stream.  The  error  comes 
from  heat  conduction  at  the  ends  of  the  probe,  which  lie  in  the  nozzle 
boundary  layer.  The  Kn  =  0.  28  data  in  Figure  25  were  corrected  for  this 

CO 
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Figure  25.  Stagnation  Line  Profiles  on  Cylinders  at  =  5.  4  in  Argon 

effect  and  both  corrected  and  uncorrected  profiles  are  shown  there.  The 
correction  makes  use  of  the  measured  point  temperatures  at  the  thermo¬ 
couple  junctions  of  the  probe  and  of  the  analytical  solutions  in  NACA  TN 
2599  by  Scadrou  and  Warshawsky  (Reference  13).  These  corrections 
require  approximate  values  of  temperature  and  velocity  in  the  shock 
wave.  These  were  obtained  from  the  solutions  of  Liepmann,  et  al  (Refer¬ 
ence  18).  The  corrected  curve  at  =  0.28  is  also  shown  in  Figure  33 
elsewhere  in  the  report  where  it  is  compared  to  a  Monte  Carlo  calculation. 

Shov/n  in  Figure  26  are  stagnation  line  profiles  at  various  Knudsen 
numbers  in  argon  and  nitrogen.  In  Figure  27  can  be  seen  the  development 
of  a  boundary  layer  between  the  shock  and  body,  as  discussed  above,  as 
the  Knudsen  number  decreases  from  0.  098  to  0.  0062. 

2.  4,  2.  2  The  Flow  about  Wedges 

As  in  the  case  of  the  cylinder,  the  flow  field  in  the  stagnation  region 
of  various  wedges  was  investigated. 
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Figure  26.  Stagnation  Line  Profiles  on  a  0.  0153  in.  Diameter 
Cylinder  in  Nitrogen  and  Argon 
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Figure  27.  Stagnation  Line  Profiles  on  Cylinders  at 
M  =  6.  3  in  Nitrogen 


In  Figure  28  is  shown  recovery  temperature  profiles  in  the  flow 
field  about  a  sharp  (nose  radius  less  than  0.  0003  inch)  20°  symmetrical 
wedge.  It  can  be  seen  that  a  disturbance  extends  about  12  mean  free 
paths  in  front  of  the  wedge.  As  one  might  expect  the  strength  of  the  dis¬ 
turbance  is  very  slight.  It  can  be  seen  in  Figure  29  how  the  disturbance 
increases  in  strength  as  the  mean  free  path  becomes  smaller. 

In  Figures  30  and  31  the  flow  field  in  front  of  a  blunt  18°  included 
angle  wedge  is  presented  for  a  warm  and  liquid-nitrogen  cooled  body. 

As  one  might  expect,  the  beginnings  of  the  disturbance  moves  in  toward 
the  body  as  it  is  cooled.  A  similar  effect  can  be  seen  in  the  curves  of 
ligure  32.  Here  we  can  see  the  entire  shock  structure  pulled  in  toward 
the  body  for  the  cold  cases,  and  the  large  drop  in  the  temperature  through 
the  boundary  layer. 


Figure  28.  Lines  of  Constant  Wire  Temperature  in  the  Stagnation  Region 
of  a  Sharp  20°  Wedge  at  =  6.7  in  Nitrogen 
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Figure  29.  Stagnation  Line  Profiles  on  a  Sharp  20°  Wedge 
at  M  ~  6.  5  in  Nitrogen 


Figure  3U.  Lines  of  Constant  Wire  Temperature  in  the  Stagnation  Region 
of  a  Blunt  18°  Wedge  at  M  =  6.7  in  Nitrogen 


2.  5  COMPARISONS  OF  MONTE  CARLO  THEORY  WITH  EXPERIMENT 


2.  5.  1  Freft  Molecule  Equilibrium  Temperature  Measurements  and  Cal¬ 
culations _ _ 

The  operation  of  the  equilibrium  temperature  probe  and  the  flow 
field  survey  performed  with  it  in  the  TRW  Low  Density  Wind  Tunnel  were 
discussed  in  the  preceding  section.  Also  mentioned  previously  was  the 
fact  that  the  Monte  Carlo  programs  can  compute  the  energy  flux  and 
molecular  number  flux  impinging  upon  a  circular  region  of  arbitrary  size 
and  location  in  the  flow  field.  Molecules  do  not  collide  with  the  circle; 
they  pass  through  it  in  their  normal  movement  in  the  field;  hence,  there 
is  no  disturbance  to  the  flow.  These  circles  simulate  equilibrium  temper¬ 
ature  probes,  and  by  assuming  a  reflection  process,  the  recovery  temper¬ 
ature  which  an  actual  probe  at  that  location  would  reach  can  be  computed. 

Assuming  complete  accommodation  and  diffuse  reflections  and 
defining: 

e  =  energy  flux  impinging  on  the  probe 
n  =  molecular  number  flux  to  the  probe 
Tw  =  probe  recovery  temperature 
k  =  Boltzmann  constant 

An  energy  balance  for  a  monatomic  gas  yields: 
e  -  2  n  kTw 

(Since  molecules  effusing  from  the  wire  in  Maxwellian  equilibrium 
carry  an  average  energy  per  molecule  of  2  kTw.  ) 

Then: 

Tw  f 

Tta  “  2  k  n  TtB 

The  Monte  Carlo  programs  also  compute  the  distribution  function 
from  which  the  probe  recovery  temperature  can  be  computed.  However, 
a  large  amount  of  numerical  data  handling  would  be  required  to  do  this, 
and  the  above  procedure  is  far  simpler.  Note  that  since  our  "ideal" 
probes  in  the  computer  program  do  not  disturb  the  flow  they  can  be  made 
larger  than  physical  probes  could  be,  the  only  requirement  being  that  they 
be  small  compared  to  gradients  in  flow  properties.  Most  of  the  calcu¬ 
lations  of  probe  temperature  which  will  be  discussed  next  were  made  with 
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probes  of  diameter  d  =  0.  6X.CQ.  At, the  point  in  the  shock  layer  where 
gradients  were  largest,  check  calculations  were  made  with  probes  of 
diameter  d  =  0.  and  0.  2\cq.  There  was  negligible  effect  of  "ideal" 
probe  size  on  the  recovery  temperature. 

In  Figure  33  a  comparison  is  shown  of  measured  and  calculated 
recovery  temperature  profile  on  the  stagnation  line  of  a  cylinder  at 
Kn  =  0.  30.  The  agreement  is  quite  good  in  general,  but  the  profiles  differ 
near  the  free  stream  and  near  the  body.  In  the  calculation,  the  body  was 
held  at  the  free  stream  total  temperature,  T^  ;  while  the  cylinder  in  the 
wind  tunnel  took  on  the  recovery  temperature,  which  at  this  Knudsen 
number  is  somewhat  higher  than  Tt  .  This  is  the  reason  for  the  drop  off 
in  computed  temperature  near  the  body. 

The  free  stream  end  of  the  computed  profile  does  appear  to  be 
rather  too  rapid  a  beginning  of  the  disturbance,  resembling  a  Navier- 
Stokes  profile,  while  the  measured  profile  resembles  H.  W.  Liepmann's 
Krook  model  calculations  which  are  thought  to  be  close  to  the  truth.  This 
matter  is  still  being  studied.  The  determination  of  the  mean  free  path 
for  argon  in  the  density-temperature  regime  used  for  the  measurements 
involved  extrapolation  of  existing  viscosity  data,  and  there  may  be  an 
error  introduced  here.  More  comparisons  of  this  kind  will  be  made  in 
tile  near  future. 

2.  5.  2  Density  Measurements  and  Calculations  on  a  Sphere 

Some  electron-beam  density  measurements  on  the  stagnation  line 
of  a  sphere  at  Kri  =  0.  03,  S  =  3.  5  were  kindly  made  available  to  us  by 
Dr.  David  Russell  of  JPL.  A  comparison  of  these  measurements  with 
Monte  Carlo  results  is  shown  in  Figure  34.  The  calculation  was  made 
on  the  7040  computer  by  Professor  Bird  in  Sydney  as  a  test  case  for 
some  future  work.  This  result  is  not  as  well  converged  as  some  of  our 
other  results,  and  there  is  some  scatter.  However,  it  is  clear  that  there 
is  very  good  agreement  in  the  density  profile  except  near  the  body  whore 
the  calculated  density  is  slightly  lower.  This  is  not  unexpected  since  the 
density  calculation  represents  a  mean  value  for  a  small  region  of  space 
slightly  off  the  stagnation  line.  Away  from  the  body,  this  is  a  good  repre¬ 
sentation  of  the  axis  value  since  the  velocity  normal  to  the  axis,  v  veloc¬ 
ity,  is  negligible  there.  Near  the  body,  the  v  velocity  increases,  however, 
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Figure  33.  Comparison  of  Measured  and  Calculated  Free  Molecule 
Equilibrium  Temperature  Along  Stagnation  Line  of  a 
Cylinder 


Figure  34.  Comparison  of  Measured  and  Calculated  Density  Along 
Stagnation  Line  of  a  Sphere 
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and  our  spatial  average  density  does  not  rise  as  much  as  the  true  stag¬ 
nation  line  density.  Reference  to  a  picture  of  typical  cell  arrangement, 
say  Figure  8,  will  make  this  clearer. 

It  is  planned  to  do  the  calculation  again  with  smaller  cells  in  the 
stagnation  region  and  a  larger  sample  size. 

2.  5.  3  Measured  and  Calculated  Drag  on  Cylinders  and  Spheres 

A  comparison  of  Monte  Carlo  drag  calculations  on  cylinders  with  a 
collection  of  measured  drag  data  from  the  published  literature  is  shown 
in  Figure  35,  covering  the  range  from  Kn  =  0.  025  to  100,  at  S  =  5.  The 
measurements  were  made  in  air  for  body  temperatures  ranging  from 
Tb/Tjo  =  8.  0  (free  stream  total  temperature)  to  T^/T^  =  9.  15,  which  is 
close  to  free  molecule  recovery  temperature  for  a  diatomic  gas„  The 
calculations  are  for  a  body  at  T^/T^  =11,  which  is  free  stream  total 
temperature  for  a  monatomic  gas  at  S  =  5.  This  temperature  difference 
should  be  negligible  at  the  lower  Knudsen  numbers,  but  it  does  have  a 
significant  effect  in  the  free  molecule  results.  The  calculated  CD  shown 
at  Kn  =  100,  Cj-j  =  2.  99,  is  within  0.  3%  of  the  correct  theoretical  value 
for  T^/Tgj  =11,  which  is  3.  0.  The  diffuse  free  molecule  limit  shown  is 
the  correct  theoretical  value  for  T^/T^  =  8.  0.  An  asymptotic  expression 
for  Cj}  free  molecule  is 


for 

S  >  5 

from  which  these  results  may  be  verified. 

The  maximum  expected  experimental  error  is  shown  in  the  figure. 
Estimates  of  the  statistical  uncertainty  in  the  Monte  Carlo  results  are 
shown  in  Table  2.  They  range  from  0.  5%  to  2.  3%. 

A  similar  comparison  for  sphere  drag  at  M  =  4  and  6  is  shown  in 
Figure  36,  where  Ren  is  used  as  parameter.  The  TRW  measurements 
of  sphere  drag  obtained  during  this  investigation  are  included  in  this 
comparison. 


-37- 


LEGEND! 


—  0DIFFUSE  FREE  MOLECULE  FLOW  LIMIT  T^T*  =  8.0,  M  =  5.92 

^^TWO-DIMENSIONAL  STRIP,  WILLIS  REF  7 

O  MASLACH,  1963  REF  8  M  =  5.92,  T^/T*  =  8.0  to  9.15,  AIR 

■  MC 

I  “ 

Lit 

3NTE 

TIMAT 
AITS  S 

CARLO 

ED  EXP 
EE  TA81 

CALCULA 

ERIMENTA 
E  2 

TIONS  M 

L  ERROR  L 

=  5.4 

IMITS 

8,  T,/T 

FOR 

00=  II,  w 

MONTE  C 

ONATOM 

ARLO  ERR 

1C  G/ 

OR 

LS 

l 

jDL-n 

3 - 

O 

0  o 

—  ' 

r~  * 

f 

/ 

5 - 

n  0 

O 

O  I 

0 

1 

3  - 

o 

)  ■  O 

>1 

oc 

1 

Q - 

■ 

1 

— 

0.02  0.04  0.06  0.1  0.2  0.4  0.6  0.1  2  4  6  10 


Figure  35.  Comparison  of  Measured  and  Computed  Cylinder  Drag 
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Figure  36.  Comparison  of  Measured  and  Computed  Sphere 
Drag 


The  agreement  between  measured  and  computed  drag  on  these  body 
shapes  over  the  range  of  Kn  and  M  considered  is  excellent. 

More  comparisons  are  planned  at  other  Knudsen  numbers,  Mach 
numbers  and  body  temperatures,  and  for  other  body  shapes  such  as 
wedges  and  cones. 

2.6  CONCLUSIONS 

A  Monte  Carlo  method  for  computing  flow  fields  and  body  forces 
has  been  extensively  developed.  Monte  Carlo  calculations  have  been 
ma.de  for  cylinders,  spheres,  sharp  wedges,  and  sharp  cones  over  the 
Knudsen  number  range  0.025  to  100.0  and  Mach  number  range  3.8  to  11.0. 
These  results  are  unique  in  that  there  are  no  other  theoretical  results 
for  complete  bodies  over  the  entire  transition  regime. 

Drag  measurements  on  spheres  and  slightly  blunted  slender  cones 
have  been  made  in  the  TRW  Low  Density  Wind  Tunnel  over  the  Knudsen 
number  range  0.  005  to  0.  025  and  Mach  number  range  4.  0  to  6.  0.  The 
agreement  with  existing  drag  data  in  the  literature  is  very  good. 

Flow  field  surveys  with  an  equilibrium  temperature  probe  have  been 
carried  out  in  the  TRW  Low  Density  Wind  Tunnel  on  cylinders  and  wedges 
in  the  Knudsen  number  range  0.  005  to  0.  30  and  Mach  number  range  4.  0 
to  6.  0.  These  results  provide  information  in  a  Knudsen  number  range 
not  previously  studied. 

A  number  of  comparisons  of  the  Monte  Carlo  theoretical  results 
with  experimental  results  has  been  made;  viz, 

1)  Measured  and  computed  equilibrium  temperature  pro¬ 
files  on  the  stagnation  line  of  a  cylinder  at  Kr  =  0.  3. 

2)  Measured  and  computed  density  profiles  on  the  stag¬ 
nation  line  of  a  sphere  at  Kn  =  0.  03. 

3)  Drag  measurements  and  calculations  on  cylinders  and 
spheres  over  range  Kn  =  0.  025  to  100. 

In  all  comparisons  above,  the  agreement  was  quite  close,  strongly 
suggesting  that  the  Monte  Carlo  theory  can  describe  the  transitional  flow 
regime  with  considerable  accuracy. 
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2.  7  RECOMMENDATIONS  FOR  FUTURE  WORK 


While  the  degree  of  agreement  to  date  between  measurements  and 
the  predictions  of  the  Monte  Carlo  method  is  very  encouraging,  there  ere 
many  more  comparisons  which  can  be  and  should  be  made;  e.  g.  , 

1)  comparisons  with  the  TRW  equilibrium  temperature 
surveys  off  the  stagnation  line  and  at  other  Knudsen 
numbers; 

2)  comparison  with  measurements  of  heat  transfer  to  fine 
wires,  a  large  collection  of  which  exists  in  the  literature; 

3)  it  has  been  claimed  that  the  distribution  function  of  a 
gas  can  be  measured  by  an  electron  beam  technique. 

Such  measurements,  if  available,  could  provide  another 
direct  check  on  the  calculated  distribution  function. 

4)  extensive  drag  measurements  exist  in  the  literature, 
especially  on  spheres  and  cones.  More  direct  compar¬ 
isons  should  be  made  of  measured  and  computed  drag. 

The  Monte  Carlo  method  presently  uses  the  rigid  elastic  sphere 
collision  model  with  a  perfectly  accommodating,  diffusely  reflecting  sur¬ 
face  model.  Other  collision  and  surface  reflection  models  should  be 
investigated. 

Extension  of  the  method  to  include  internal  degrees  of  freedom  locks 
to  be  within  the  realm  of  possibility  and  this  should  be  pursued. 

The  method  can  very  easily  be  extended  to  include  other  body 
shapes;  e.  g.  ,  wedges  and  cones  of  arbitrary  bluntness  and  apex  angle. 

Comparison  with  experimental  data  for  a  very  complicated  flow 
which  can  be  handled  inexpensively  on  the  computer  should  be  made.  An 
experimental  program  at  California  Institute  of  Technology  to  measure 
shock  reflection  from  a  wall  is  underway.  This  flow  could  be  described 
inexpensively  on  the  computer  since  it  is  one-dimensional,  and  it  would 
be  a  very  exacting  test  of  the  theory. 
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SECTION  III 


WAKE  TURBULENCE  MEASUREMENTS 

(HYPERSONIC  WAKE  AND  BOUNDARY  LAYER  TECHNOLOGY) 

3.1  INTRODUCTION  AND  BACKGROUND 

A  knowledge  of  the  fluid  dv.^amical  turbulent  wake  properties  is 
needed  to  understand  and  predict  radar  scattering  from  reentry  vehicle 
wakes.  Among  the  turbulent  wake  properties,  the  three-dimensional 
correlation  length  plays  a  major  role  in  the  current  theories  of  radar 
backscatter.  Previous  experiments  in  the  TRW  Systems  ballistic  range 
measured  the  wake  turbulence  behind  .22  caliber  proj  ictiles  traveling 
at  4000  feet  per  second.  A  single  anemometer  probe  was  placed  close 
to  the  projectile  path  to  sense  the  quickly  growing  wake  behind  the  pro¬ 
jectile;  analysis  of  the  fluctuating  data  providing  correlation  times. 
Radars,  however,  sense  space  correlation  lengths.  Therefore,  it  was 
an  objective  to  obtain  experimental  data  of  wake  turbulence  fluctuations 
at  more  than  one  point  simultaneously.  It  was  an  additional  objective 
to  relate  the  measured  space  correlation  lengths  to  those  in  the  turbulent 
wake  behind  a  typical  slender  reentry  vehicle.  The  accomplishments 
during  the  contract  period  are  described  in  the  following  sections. 

3.  2  SPACE  CORRELATION  LENGTH  MEASUREMENTS 

This  program  was  begun  with  two  hot-film  anemometers  placed 
adjacent  to  the  projectile  path  in  a  newly  established  range.  Two  types 
of  projectiles  were  launched  from  a  .22  caliber  rifle  at  4000  fps,  one 
type  being  more  slender  than  the  other.  The  fluctuating  anemometer 
data  was  produced  shortly  behind  the  projectile  by  the  growing  wake 
which  enveloped  the  anemometer.  Later  statistical  analysis  of  the  re¬ 
corded  data  on  a  digital  computer  produced  auto  and  cross  correlations 
together  with  their  Fourier  transforms,  power  spectral  densities  and 
cross  spectra.  Eddy  convection  velocities,  space  correlation  length 
scales,  and  time  scales  were  also  produced,  the  latter  having  been 


*The  sphere-ogive -cylinder  shapes  of  these  commercial  bullets  is 
shown  in  Figure  4  of  the  report  identified  at  the  end  of  Section  3.  2. 
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demonstrated  earlier  in  a  measurement  program  with  a  single  anemo¬ 
meter  (19).  Space  correlation  length  of  the  wake  turbulence  was  shown 
for  the  first  time  in  these  measurements  to  increase  with  distance  be¬ 
hind  the  projectile  as  predicted  by  the  Lees-Hromas -Webb  (20,  21) 
theory.  Previous  measurements  of  the  wake  core  turbulence  (22,23) 
had  shown  a  constant  space  correlation  length  up  to  9000  diameters 
behind  projectiles.  These  latter  measurements  had  been  made  using 
optical  methods  (Schlieren  and  interferometer)  which  integrate  along 
the  light  beam  passing  completely  through  tue  wake.  In  contrast,  the 
anemometer  makes  a  local  measurement,  its  length  (.  20-inch)  being 
much  smaller  than  the  wake. 

These  measurements  were  presented  at  the  AIAA  Fifth  Aerospace 
Sciences  Meeting,  New  York,  January,  1967  as  AIAA  Paper  No.  67-19. 
They  are  also  available  as  TRW  Report  No.  07854-6021-R000,  BSD- 
TR-67-82.  This  report  has  been  scheduled  for  publication  in  the 
February,  1968,  issue  of  the  AIAA  Journal. 

1 

3.  3  EXTENDED  X/D  DATA 

The  next  series  of  measurements  utilized  the  more  blunt  .22  caliber 
projectile  and  two  different  spacings  between  the  anemometers.  A 
two- data  estimate  was  thus  obtained  for  the  space  correlation  function 
and  its  integral,  the  space  correlation  length.  Three  different  anemo¬ 
meter  temperatures,  62°C,  166°C  and  270°C,  were  used  to  provide  data 
for  reduction  to  velocity,  density,  and  temperature  fluctuations.  In  all, 
data  from  35  shots  were  reduced.  The  statistical  analyses  were  made 
over  a  much  greater  range  of  X/D  than  that  used  previously.  As  shown 
in  Figure  37,  the  data  from  X/ D  =  100  to  160,000  was  used  in  the  analyses 
The  delay  times  in  the  auto  and  cross  correlations  of  Figure  37,  as  well 
as  the  frequencies  in  power  spectral  densities  of  Figure  38,  extend  over 
four  orders  of  magnitude.  The  definitions  of  autocorrelation  and  cross 
correlation  coefficients  used  in  Figure  37  were  retained  from  the  pre¬ 
viously  mentioned  report,  namely, 


C  (At)  = 

g 


*  Bluntness  ratio  of  0.  6. 


+  At) 
g?t] 
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POWtR  SPECTRAL  DENSITY  (AM  UNITS) 


and 


C  ,  (Ax,  At)  =  - Z_£.(°»t)  h(Ax,t  +  At) _ _ 

8  |z  g(t)  g(t)J  A/  2  h(Ax,  t)  h(Ax,  t)J  11  2 


where  g(t)  and  h(Ax,  t)  are  the  digitized  data  from  the  two  anemometers, 
the  respective  mean  values  having  been  previously  subtracted  away. 

^he  time  integral  scale,  T,  which  is  shown  in  Figure  39  factored  by 
the  projectile  velocity  U^produced  an  acceptable  trend  over  the  length 
of  the  wake.  The  statistical  techniques  used  to  obtain  these  quantities 
are  described  in  the  previously  mentioned  reports. 

3.4  OFF-AXIS  DATA 

Five  of  the  previously  mentioned  35  shots  that  were  reduced  utilized 
a  near-axis  location  on  one  anemometer  and  an  off-axis  location  on  the 
other.  This  configuration  was  tested  in  response  to  a  SAMSO/ Aerospace 
request  for  an  off-axis  survey.  The  planning  for  the  off-axis  shots  is 
snown  in  Figure  40  where  the  wake  edge  and  anemometer  probe  locations 
are  shown  in  units  of  projectile  radii  R.  At  Y/R  =  1,  the  probe  and 
projectile  make  contact.  The  off-axis  probe  was  one  projectile  diameter 
farther  away  from  the  projectile  axis  than  the  near-axis  probe. 

The  redus„ed  data  that  was  obtained  was  autocorrelation,  power 
spectra,  and  time  integral  scales;  probe  separation  being  too  great  for 
space  correlation  data.  The  off-axis  time -correlation  data  followed  the 
same  trends  as  the  near-axis  data,  which  is  shown  in  Figures  37,  38,  39. 
No  additional  off-axis  tests  were  conducted  during  the  contract  period. 

The  "theory"  wake  edge  shown  in  Figure  40  was  provided  by  a  new 
calculation  using  the  Lees-Hromas-Webb  theory  together  with  a 
CD  =  °*4,  which  is  more  accurate  for  this  projectile  than  the  previously 
used  CD  =  0.7.  (This  revision  is  reflected  in  all  of  the  theoretical 
results  shown  in  this  section.)  In  the  theory,  the  wake  width  is  really 
an  energy  diffusion  width,  whereas  the  experimental  wake  width  is  the 
average  across  the  scalloped  wake  projection  on  a  shadowgraph.  In 
view  of  these  limitations  of  the  comparison,  the  agreement  in  Figure  40 
between  the  Lees-Hromas-Webb  theory  prediction  and  the  measured 
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wake  edge  is  considered  satisfactory. 

3.  5  TWO-DATA  SPACE  CORRELATION  LENGTH  SCALE 

Anemometer  probe  separations  Ax  of  R/  2  and  R  were  used  in  ac¬ 
quiring  data.  Cross  correlation  output  produced  the  statistical  data 
needed  to  estimate  the  space  correlation  function,  the  integral  of  which 
is  the  space  correlation  length  (distance  over  which  fluctuations  are  cor¬ 
related).  The  cross  correlation  data  were  first  averaged  over  10  to  15 
shots  with  all  probe  temperatures  being  utilized.  The  averaged  cross 
correlations  were  then  smoothed  in  X/ D  with  a  weighted  average  scheme; 
50  percent  of  the  point  in  question  was  averaged  with  25  percent  of  each 
of  two  adjacent  data  points.  A  curve  was  fit  to  the  two  data  for  Ax  =  R/  2 
and  R  at  each  X/D.  This  curve  extended  from  unity  at  a  zero  separation 
(Ax  =  0)  to  zero  at  a  large  separation  (Ax  =  °°).  The  integral  of  that 
curve,  namely,  the  space  correlation  length  A,  is  shown  in  Figure  41. 

It  agrees  with  the  prediction  of  the  Lees -Hromas -Webb  theory,  thus 
verifying  one  of  the  major  features  of  the  theory  that  affect  the  radar 
backscatter  cross  section. 


AS- WEBB  THEORY 

^“NEW  TUyy 
-  OLD  & 

100  1000  10,000  100,000  1,000,000 

X/D 

Figure  41.  Space  Correlation  Length 
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oeogr ate  studies  of  A  measurements,  considering  the  low  anemometer 
temperature  data  and  high  anemometer  temperature  data  separately, 
gave  no  trend  of  A  with  probe  temperature.  The  random  irregularities 
of  the  A  data  were  greater  in  this  case  because  a  smaller  number  of 
shots  contributed  to  the  averaging.  Thus  any  conclusions  about  the  lack 
of  trend  in  the  data  must  be  qualified  by  the  averaging  irregularities.  It 
will  be  explained  in  a  later  section  that  the  data  taken  at  low  anemometer 
temperature  are  mostly  wake  temperature  data  and  those  at  high  ane¬ 
mometer  temperature  are  mostly  wake  velocity  data.  Thus  the  lack  of 
trend  of  A  with  anemometer  temperature  indicates  that  a  single  A 
characterizes  both  temperature  and  velocity  fluctuations,  within  the 
precision  of  the  measurements. 

By  Taylor  s  hypothesis,  the  product  of  the  time  integral  scale  T 
and  the  mean  wake  velocity,  as  predicted  by  the  Lees-Hromas-Webb 
theory  is  comparable  to  A  .  This  comparison  is  indeed  shown  to  be  good 
for  both  the  old  data,  which  was  obtained  during  the  first  program  des¬ 
cribed  in  3.2,  and  the  new  extended-range  data.  The  lesser  X/ D  range 
obtained  for  A/D  compared  with  TUw/D  suggests  that  extraneous  effects 
in  the  ballistic  range  alter  space  correlations  more  quickly  than  they 
alter  time  correlations. 

3.6  MEAN  VELOCITY,  TEMPERATURE  AND  DENSITY 

The  original  reduction  of  mean  wake  velocity  and  temperature  was 
taken  from  oscilloscope  photographs,  this  reduction  being  later  confirmed 
by  digital  computer  output.  The  data  were  smoothed  and  averaged  between 
three  to  five  shots.  Data  for  two  different  probe  temperatures  were  thus 
obtained  fromprojectile  wakes  and  compared  with  calibrated  responses 
obtained  with  the  probes  in  a  heated  jet  and  in  an  oven.  The  data  obtained 
with  low  probe  tempersture  was  mostly  temperature,  whereas  the  high 
probe  temperature  data  was  mostly  velocity. 

Mean  wake  velocity  and  temperature  obtained  by  this  technique  is 
shown  in  Figures  42  and  43.  The  density  can  be  obtained  directly  from 
the  temperature  data  because  the  pressure  is  constant  in  the  far  wake, 
and  perfect  gas  laws  are  valid.  Both  the  velocity  and  temperature  are 
somewhat  lower  than  the  Lees-Hromas-Webb  theory  prediction. 
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As  shown  in  Figure  40,  the  near-axis  data  are  located  a  significant 
fraction  of  the  wake  radius  away  from  the  wake  axis  when  X/ D  is  small. 
In  order  to  make  a  better  comparison  with  the  axis  values  that  are  pre¬ 
dicted  by  the  Lees-Hromas-Webb  theory,  an  equivalent  to  the  measured 
data  was  constructed  using  a  profile  curve  fitted  to  the  wake  width  which 
passed  through  the  near-axis  data.  The  axis  values  of  the  fitted  curves 
are  shown  in  Figures  42  and  43  as  the  axis  equivalents  of  the  measured 
data.  The  axis  velocity  gives  good  agreement  with  the  Lees-Hromas- 
Webb  prediction  up  to  the  limit  of  calibration  at  X/D  =  2,000.  This  agree 
ment  verifies  another  major  feature  of  the  Lees-Hromas-Webb  theory 
which  is  directly  related  to  the  radar  Doppler  return.  Somewhat  less 

agreement  is  obtained  by  the  axis  temperature  in  Figure  43,  although  the 
correct  trend  is  obtained. 


Greater  X/D  range  was  obtained  with  mean  convection  velocities  U  . 
that  is.  the  velocity  of  convection  of  fluctuations  from  one  anemometer  to 
the  other,  as  determined  by  best  cross  correlation  delay  time.  The  av¬ 
erage  of  at  least  9  and  as  many  as  25  shots  are  shown  in  Figure  44.  It 
is  reasonable  that  the  convection  velocity  Uc  falls  below  th,  previously 
determined  mean  wake  velocity,  shown  in  Figure  42,  bee,  »  the  average 
over  many  more  shots  was  used  for  the  convection  velocity.  The  probe- 
to-axis  distance  was  greater  for  most  shots  than  those  few  use  to  obtain 

the  mean  wake  velocity  measurement,  resulting  in  a  higher  mea,  wake 
velocity  than  convection  velocity  up  to  the  limit  of  calibration.  I,  s 
indeed,  reassuring  that  a  trend  parallel  to  the  predicted  velocity  was 
obtained  for  convection  velocity  up  to  an  extended  X/D  of  70,  000. 


3.  7 


FLUC  TUATION  MEASUREMENTS 
TEMPERATURE  AND  DENSITY 


VELOCITY, 


Among  the  35  shots,  data  was  taken  with  three  different  probe 
temperatures  so  that  the  fluctuations  of  velocity,  tempo  ,Ure  and 
density  could  be  obtained.  The  calibration  function  q/'l  *.  +  Bx/U  was 
utilized  with  mean  and  fluctuating  components  being  substituted  and 
then  mean  squared.  The  resulting  expression  contained^  and  the 
mean  properties  as  known  factors,  and  WZ,  U'T',  and  T^as  unknowns 
The  three  sets  of  data  formulated  three  linear  equations  that  were  in¬ 
verted  to  obtain  the  three  unknowns.  A  qualified  success  was  obtained 
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Figure  43.  Mean  Wake  Temperature 


CONVECllGN  VELOCITY 


Figure  44.  Convection  Velocity  Compared  with  Mean  Wake  Velocity 


at  X/ D  =  840  where  the  following  values  were  obtained 
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where  the  density  fluctuations  were  obtained  from  the  temperature  fluc¬ 
tuations  in  much  the  same  manner  as  that  used  for  the  mean  values. 


These  values  are  in  remarkable  agreement  with  Townsend's  (24) 
cylinder  results  and  the  recent  axi- symmetric  measurements  of  Demetri- 
ades  (25)  and  Gibson  (26).  The  necessary  reduction  of  the  cross  product 
produced 


TFT7  _  ,  14 

1/2  — -TTT  ~  ~1’  16 
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which  has  the  correct  sign  and  is  only  slightly  too  large.  At  greater  X/D 
the  agreement  in  magnitude  was  poorer,  although  the  correct  signs  were 
maintained.  These  results  definitely  prove  the  feasibility  of  the  measure¬ 
ment  technique  and  serve  to  indicate  the  needed  refinements.  These  re¬ 
finements  include  a  wider  range  of  probe  operating  temperatures,  up  to 
300°C,  and  calibrations  to  lower  velocities,  down  to  15  fps. 

The  underlying  philosophy  of  the  Lee s-Hromas -Webb  theory  con¬ 
siders  the  turbulence  structure  of  the  far  wake  as  basically  similar  to 
low  speed  turbulent  structure.  The  verification  obtained  in  these  mea¬ 
surements  supports  that  philosophy  and  lends  credence  to  its  use  in  des¬ 
cribing  electromagnetic  backs catter.  To  the  extent  that  the  electron  field 
can  be  considered  locked  to  the  flow  field  of  the  R/ V  wake,  the  understanding 
of  electromagnetic  backscatter  has  been  enhanced,  especially  with  respect 
to  space  correlation  length  A  and  the  wake  velocity.  These  quantities  con¬ 
tribute  to  the  radar  cross  section  and  Doppler,  respectively.  Both  quan¬ 
tities  have  been  shown  to  be  predictable  by  the  Lees-Hromas-Webb  theory. 
Thus,  R/V  wake  predictions  using  this  theory  have  support  from  controlled 
laboratory  measurements. 

3.8  SPHERE  WAKE  MEASUREMENTS 

Half-inch  spheres  were  launched  with  near-zero  scatter  at  6500  fps, 
using  a  new  smooth-bore  gun.  A  wake  measurement  program  was  begun 
to  better  resolve  velocity  and  temperature  fluctuations.  Higher  probe 
temperatures  (to  450°C),  a  greater  range  of  probe  temperatures  (300°C), 
and  a  refined  jet  velocity  instrumentation  for  calibration  were  employed 
in  this  program.  An  interruption  in  the  program  was  allowed  for  the  cone 
wake  measurements,  and  the  program  was  again  underway  at  the  end  of 
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the  contract  period. 

3.  9  CONE  WAKE  MEASUREMENTS 

Over  100  launchings  of  a  12-1/2  degree  half-angle  cone  with  3/ 8-inch 
base  diameter  were  accomplished,  using  a  rifled  barrel  to  obtain  speeds 
of  6500  fps.  Measurements  were  obtained  using  probe  temperatures  of 
450°C.  A  range  pressure  of  100  mm  was  used  to  permit  sabot  separation 
with  minimum  scatter  of  the  projectiles.  Data  records  from  eight  laun¬ 
chings  were  delivered  to  the  TRW  Computer  and  Data  Reduction  Center. 
The  reduction  and  statistical  analysis  of  the  data  was  underway  at  the  end 
of  the  contract  period.  From  these  data,  the  correlation  length  scales 
will  be  obtained  that  are  comparable  with  a  slender  reentry  vehicle. 

3.10  CONCLUSIONS 

1.  The  understanding  <f  electromagnetic  scattering  from 
reentry  vehicle  wakes  has  been  enhanced  by  the  verif¬ 
ication  of  the  Lees-Hromas-Webb  theory  under  labora¬ 
tory  conditions.  The  measurements  of  space  correlation 
length  and  wake  velocity  were  adequately  predicted  by 
the  theory  these  quantities  being  especially  significant 
for  radar  backs catter. 

2.  Space  correlation  length  for  temperature  fluctuations 
is  the  same  as  that  for  velocity  fluctuations  within 
the  precision  of  the  data. 

3.  Time  correlation  measurements  can  be  made  to  much 
greater  X/D  than  that  for  space  correlation  measurements. 

4.  Time  correlation  scales  show  no  effect  of  off-axis  location. 

5.  Comparison  of  the  results  to  those  in  the  turbulent  wake 
behind  a  slender  reentry  vehicle  will  be  possible  when 
che  data  from  the  slender- cone  program  is  reduced. 

6.  Mean  wake  velocities,  densities  and  convection  velocities 
»n  the  wake  can  be  measured  with  relative  ease  using  ane¬ 
mometers  in  the  ballistic  range.  Fluctuating  components 
of  velocity  and  density  are  more  difficult  to  obtain. 

3.11  RECOMMENDATIONS 

1.  The  measurements  of  space  correlation  length  and  velocity 
and  mass  density  fluctuations  using  anemometers  should 
be  extended  to  wakes  behind  projectiles  with  higher  speeds, 
higher  Mach  number,  larger  diameter  and  varying  shape. 

More  complete  comparison  with  slender-cone  reentry 
vehicle  wakes  requires  the  variation  of  these  important 
parameters. 
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?,.  Electron  density  fluctuations  should  be  measured  at  the 
same  time  as  the  mass  density  fluctuations.  Radar  back- 
scatter  measurements  from  reentry  vehicle  measurements 
can  then  be  compared  with  the  mass  density  measurements 
and  electron  density  measurements  available  from  the 
ballistic  range. 
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SECTION  IV 

4.  I  INTRODUCTION  AND  BACKGROUND 

The  objective  „£  this  task  i.  to  develop  an  analytical  description  of 
the  compressible,  turbulent,  chemically  reacting  boundary  layer  and  to 
confirm  the  essential  features  of  the  theoretical  model  by  experimenta¬ 
tion.  The  primary  reason  for  performing  this  task  is  to  provide  the  nec¬ 
essary  initial  conditions  for  the  calculation  of  the  near  wake.  These 
calculations  require  definition  of  the  dynamic,  thermodynamic  and  chem- 
ical  state  of  the  boundary  layer  leaving  the  body. 

The  technical  approach  to  this  formidable  problem  has  been  to: 

*  fu:w:p„t\‘tr;ciaIyde:rription  °f  an  <*•»»«>«■  ■«■!. 

*  formation' to’S'clTd10”  Coles'  mathematical  trans- 

include  pressure  gradient,  heat  transfer 

s  transfer  in  order  to  account  for  the  effects 
of  compressibility;  euects 

*  ZST'tfZXr  °£  turbuUnce  °”  chemical  reaction 
rates.  (It  has  been  assumed  that  the  energy  released 

f4d  m^cS  ""T'l0113  18  suff^iently  smallthat  the 

a«ect?d.  ,  a‘  S‘ate  °f  the  b0Undary  ^er  will  not  be 

Emphasis  has  been  placed  on  developing  a  sound  theoretical  model  and 
not  on  quick  programming  to  obtain  meaningless  numbers. 

4.  2  LOW  SPEED  TURBULENT  BOUNDARY  LAYER 

In  order  to  gain  insight  into  the  important  mechanisms  of  turbulent 
momentum  transport,  the  case  of  a  parallel  turbulent  shear  flow  has 
been  considered  prior  to  attacking  the  more  difficult  boundary  layer  flow. 

A  self-contained  analysis  for  parallel  turbulent  flows  in  dynamic 
eqm  ibrium  has  been  derived.  The  theory  predicts  a  priori  the  existence 
the  laminar  sublayer  and.  when  account  of  viscous  dissipation  effects 
near  the  wall  is  taken,  agrees  well  with  detailed  experimental  data  for 
both  mean  and  fluctuating  quantities  without  the  need  for  ad  hoc  fitting 
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constants.  The  detailed  analysis  and  comparisons  to  experimental  data 
are  given  in  TRW  Report  07854-6024-R000  entitled  "A  New  Mixing  Length 
Theory  for  Turbulent  Flows  in  Dynamic  Equilibrium"  by  W.  H.  Webb. 

The  advantage  of  the  method  developed,  namely,  the  elimination  of 
unnecessary  empirical  constants  which  have  in  the  past  been  required  to 
describe  such  flows,  is  clear.  Examination  of  details  of  experiments  per¬ 
formed  in  the  past  shows  that  these  "constants"  are  not  actually  constant 
at,  all  but  rather  may  vary  with  Reynolds  number,  for  example  (see  paper 
by  J.  O.  Hinze  "Turbulent  Pipe  Flow"  in  The  Mechanics  of  Turbulence, 
Gordon  and  Breach,  New  York,  1964).  Thus,  an  obvious  advantage  oJ 
the  analysis  developed  is  that  variations  of  this  kind  are  accounted  for. 
Further,  since  the  model  may  be  applicable  to  boundary  layers  (as  well 
as  confined  flows  for  which  it  was  developed)  questions  which  have  arisen 
in  the  past  such  as  whether  the  empirical  "constants"  are  the  same  for 
pipe  flow  and  boundary  layer  flow  need  no  longer  concern  us.  In  addition, 
the  applicability  of. the  method  to  heat  transfer  and  mass  transfer  calcula¬ 
tions  (Deissi.,-r  has  shown  how  this  may  be  done)  may  permit  such  results 
to  be  obtained  without  further  empiricism. 

The  approach  for  the  extension  of  the  above  theory  to  the  boundary 
layer  has  been  defined.  The  physical  situation  is  believed  to  be  well 
represented  by  two  regions:  an  inner  region  in  which  the  flow  behaves 
like  a  parallel  flow,  and  an  outer  region  in  which  the  velocity  does  not 
deviate  appreciably  from  the  external  velocity,  thereby  allowing  an  Oseen 
type  approximation. 

As  in  the  model  of  Townsend,  in  the  outer  region  energy  is  extracted 
from  the  mean  flow  by  the  turbulent  shear  stress  u'v'.  The  energy  is 
converted  to  turbulent  energy  and  dissipated  mainly  in  the  inner,  nearly 
constant  stress  region.  (See  Figure  45. )  Both  regions  are  bounded  by 
layers  which  are  laminar  in  the  sense  that  laminar  shear  dominates 
turbulent  shear  in  these  layers.  The  viscous  sublayer  plays  the  role  of 
transferring  stress  to  the  wall;  the  viscous  superlayer  plays  the  role  of 
converting  nonturbulent  fluid  to  turbulent  fluid  across  an  irregular  (inter¬ 
mittent)  boundary. 

Mathematically  these  regions  are  conveniently  (albeit  approximately) 
described  by  a  Couette  flow  approximation  (constant  shear)  in  the  inner 
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Figure  45.  High  Reynolds  Number  Turbulent  Bounda 
Layer 


region  and  an  Oseen  flow  approximation  (constant  mean  velocity)  in  the 
outer  region. 


If  the  scaling  velocity  for  the  turbulent  shear  is  taken  to  be 

d 

AU  =  U(d)  f  Udy 
o 

where  d  is  the  inner  layer  height,  then,  as  for  the  case  of  equilibrium 
channel  flow,  by  representing  the  turbulent  shear  stress  by 


t.  =  f(x,  y) 


mi 


where  f(x,  y)  is  an  unspecified  function  of  x  and  y,  the  assumption  of 
constant  shear  leads  to  a  solution  for  f(x,  y)  which  results  in  the  nor¬ 
malized)  shear  stress 


,  -  -(■  4) 


T(  =  o 


6b  S  y  S  d 


0  £  y  £  6 


where  6  is  the  sublayer  thickness  which  is  given  by 


6  = 
s 


VR6  6  R 

747172  6  ■  Rs  * 


<cfV 


U<6s>6s 


The  velocity  then  can  be  obtained  as 


^=Rs^[1+,°gVl/2 


where  U  is  the  friction  velocity.  This  law  is,  of  course,  of  the  form 

u  T  yuT 

jj""  =  A  +  B  log  — ~  where  the  constants  A,  B  are  related  to  the  sublayer 

Reynolds  number  Rg  .  As  shown  previously,  this  number  appears  from 
Laufer's  data  to  be  a  universal  constant. 

While  the  inner  region  is  known  to  be  similar  in  the  vs.  ™  T 

U  T  ^ 

coordinate,  the  outer  region  is  weakly  non-similar  (see  e.g.  Clauser). 

If  the  same  scaling  velocity  method  ie  also  applicable  in  this  non¬ 
equilibrium  outer  layer,  an  analytical  representation  for  this  layer  may 
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also  be  obtained.  We  again  take 

Tt  =  £(x’y>  (ir) 

for  the  outer  region.  It  is  convenient  to  let  1  -  U  =  a(x)g(q)  where 
ti  =  y/6  ,  in  which  case  the  boundary  layer  equation  gives  for  h,  where 
h'  =  g, 


h'"  |2fh"  - 


,R6)+  £ti 


h"2  =  6'  (1  -  a'6)  hh"  -  (1  -  ah')2 

a 


a  third  order  equation  in  h.  At  this  point,  we  use  the  Oseen  approxima¬ 
tion  for  the  outer  layer  so  that  h  q  and  h'  —  0,  then  calling  G  =  g'  we 
get 

f  G2  -  6'  (1  -  a'6)  qG  +  ~ 

G.  +  _D - - - 2_=  0 

na-A 

By  assuming  qua  si- similarity,  i.  e. ,  that  a'  may  be  neglected;  and  also 
by  ignoring  the  superlayer,  i.  e. ,  dropping  the  (aR^)'1  term  compared 
to  the  2fG  term,  we  get 

G.  xin  G  _ii  JL  =  o 

U  +  2f  U  2  f  U 

This  equation  may  be  solved  for  the  velocity  gradient  G  and, 
employing  the  velocity  difference 


AU  =  1 


-  /  U  dy 


as  the  scaling  velocity,  we  obtain  for  the  turbulent  shear  stress 

' 

_  ,  Vi  ^  -i. 

T*  £  ^  o  1  -  ZICjR,)1/2 

where  q  =  d/6  is  the  normalized  inner  layer  height. 

In  order  to  obtain  q  we  conserve  energy  flux  in  the  boundary  layer 
by  satisfying  the  integral  mechenical  energy  equation 
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2 


6'61  =2  /  r  ^  dn 

O 

1 

02  =  f  U(1  -  u2)  dri 


The  result  is 


t|  q2  -  (1  -  2a)  exp  (1  -  a_i) 


‘  3a 


=  1 


where  a  =  (C,R  )^2. 
'  f  s 


This  equation  has  the  solution  shown  in  Figure  46  which  indicates 
that  for  most  turbulent  boundary  layers  of  interest  where  a  *  .  1  then 
*  l  .  That  is,  the  inner  layer  is  about  10%  of  the  total  thickness  of 
the  boundary  layer,  a  result  in  accord  with  experiment. 

The  boundary  layer  thickness  is  obtained  by  satisfying  the  integral 
momentum  equation 

6*0  =  Cf 

1 

0  =  J  U(1  -  U)  dq 

o 

The  result  is  the  parametric  set 
R 

ir  =  ir  ((1  - 3b  +  bi)  exP <b  -  V  -  (*  -  3bo +  bo2)) 

Ir  =  T-  exP  <b  -  b0>  •  ct  *  (E,b2) 

o  o 

.  R6o 

bo  exp  bc  =  e  -g- 

8 

These  results  are  shown  in  Figures  47  and  48, 
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Figure  46.  Solution  for  Matching  Point,  T]#vs.  (C^R  ) 
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The  procedure  outlined  above  thus  provides  a  general  method  for 
obtaining  detailed  properties  of  the  turbulent  boundary  layer  with  a  mini 
mum  of  empirical  data.  In  order  to  obtain  quantitative  accuracy  it  is 
expected  that,  as  in  the  channel  flow  case,  the  effect  of  viscous  dissipa¬ 
tion  on  the  turbulent  shear  stress  near  the  wall  must  be  taken  into 
account.  Further,  the  effect  of  intermittency  on  reducing  the  shear 
stress  of  the  outer  layer  must  also  be  included.  These  effects  can  be 
added  as  corrections  in  a  straightforward  way  and  work  on  this  is  in 
progress. 


Of  course,  the  accuracy  of  this  method  must  be  tested  against 
experiment  sz  for  the  channel  flow  case;  however,  it  is  clear  that  if  an 
accurate  representation  of  the  turbulent  shear  stress  does  result,  then 
the  method  may  be  applied  directly  (as  shown,  e.g. ,  by  Deissler)  to  the 
calculation  of  the  transport  of  mass  or  energy  in  the  turbulent  boundary 
layer  without  further  empirical  data.  The  extension  of  this  method  to 
such  cases  as  well  as  to  cases  with  pressure  gradient,  blowing,  and 
compressibility  effects  then  holds  the  promise  of  reducing  much  of  the 
empiricism  currently  required  for  the  analysis  of  such  problems. 

4.  3  COMPRESSIBILITY  TRANSFORMATION 


The  application  of  Coles'  transformation  to  boundary  layers  with 
pressure  gradient,  heat  transfer  and  mass  transfer  has  been  accomplished 
to  the  extent  of  deriving  the  complete  set  of  Lws  of  corresponding  stations 
which  govern  these  parameters: 


cfRe  =  CfR^ 
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(previously  derived  by  Coles) 
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The  details  of  the  analysis  can  be  found  in  TRW  Report  07854-6065 -R000 
entitled  "The  Compressible  Boundary  Layer  with  Pressure  Gradient, 

Heat  Transfer  and  Mass  Transfer  and  Its  Low  Speed  Equivalent,"  by 
J.  E.  Lewis. 

4.  4  CHEMICAL  KINETICS 

4.4.  1  Hypersonic  Boundary  Layer  Kinetics 

The  first  effort  undertaken  in  the  study  of  turbulent  fluctuation 
effects  on  chemical  kinetics  was  the  determination  of  certain  laminar  or 
molecular  chemical  relations  needed  for  a  practical  application  to  the 
boundary  layer  of  an  ablating  R/V,  setting  aside  for  the  moment  the  prob¬ 
lem  of  applying  these  relations  to  turbulent  flow.  The  particular  case  of 
an  ablating  graphite  boundary  layer  was  chosen  for  its  relative  simplicity. 

The  thermochemical  model  of  dissociated  air  flowing  over  a  pure 
carbon  surface  was  selected  for  detailed  study  since,  at  surface  temper¬ 
atures  below  3000°K,  there  are  only  nine  important  electron -producing 
and  thermally  signficant  species,  (namely,  N^,  N,  O^,  O,  NO,  NO+,  e“, 
CO  and  CO^,),  gaseous  carbon  being  noteworthy  by  its  absence  at  these 
temperatures  but  rot  at  higher  temperatures. 

A  set  of  ten  homogeneous  reactions  and  three  heterogeneous  reac¬ 
tions  were  identified  through  a  literature  survey  as  the  significant 
reactions  involving  the  above  nine  species.  These  reactions  and  their 
corresponding  rate  constants  presented  in  Table  HI.  The  rates  for 
reactions  I  through  VH  were  obtained  from  Reference  28.  The  rates  for 
the  homogeneous  reactions  VHI,  IX  and  X  were  taken  from  References  29, 
30  and  31,  respectively.  The  rates  for  the  hetereogeneous  reactions  XI, 
XH  and  XIH  were  extracted  from  References  32,  33  and  34,  respectively. 

A  method  was  needed  to  transpose  these  molecular  relations  into  a 
calculation  method  to  describe  the  time  average  of  the  fluctuating  turbu¬ 
lent  flow  and  chemistry  in  the  boundary  layer.  It  was  clear  from  the 
outset  that  the  time-average  properties  of  mass -fraction,  density,  and 
temperature  could  not  be  substituted  into  the  laminar  or  molecular 
expressions  for  chemical  kinetics  to  produce  an  "effective"  turbulent 
reaction  rate.  The  effect  of  mixing,  for  example,  cannot  be  described 
by  such  a  method. 
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An  examination  of  the  Rubel-Zeiberg  (Reference  35)  theory  for  the 
effective  turbulent  reaction  rate  was  undertaken.  The  essence  of  the 
theory  is  the  substitution  cf  mean-plus -fluctuating  properties  into  the 
molect  .ar  or  laminar  kinetics  relation  and  then  taking  the  time  average 
of  the  resulting  expression.  The  new  terms  that  are  generated  are 
average  products  of  two  or  three  fluctuating  quantities  such  as  the 
density,  p;  mass  fraction  of  each  species,  a ^  ;  and  temperature  rate 
factor  of  each  reaction,  k..  .  It  is  evident  that  the  number  of  these 
average  products,  or  correlation  coefficients,  that  appear  increase 
extremely  rapidly  with  the  number  of  reactants  present  and  with  the 
number  of  important  reactions.  In  fact,  the  application  of  the  Rubel- 
Zeiberg  method  to  the  ten  homogeneous  reactions  listed  in  Table  I  would 
require  values  for  approximately  200  correlation  coefficients.  The 
magnitude  of  these  coefficien*s  could  only  be  estimated  for  a  few  of  the 
air  species  reactions  and  were  unknown  for  the  remaining  reactions. 
Efforts  to  apply  the  method  to  the  R/V  boundary  layer  ceased  and  atten¬ 
tion  was  turned  to  experimental  methods  to  obtain  the  effect  of  s.  fluctu¬ 
ating  field  on  the  chemical  kinetics.  These  efforts  included  the  definition 
of  an  experiment  to  be  conducted  in  a  wind  tunnel  and  the  studies  of  water - 
tunnel  experiments  and  suitable  low- speed  theories. 

4.  4.  2  Water-Tunnel  Kinetics  Experiments 

A  literature  search  yielded  a  series  of  three  experiments  in  a 
water  tunnel  (References  36,  38),  the  last  of  which  measured  the  effective 
rate  of  reaction  in  the  wake  of  a  cylinder.  One  opaque  reactant  was 
injected  through  the  porous  wall  of  a  cylinder,  the  other  transparent 
reactant  being  in  the  free  stream.  Reaction  rate  was  measured  by  time 
photographs  of  the  disappearing  wake.  The  calibrated  photographs  were 
measured  with  a  densitometer.  Profiles  were  thus  obtained  as  well  as 
integrated  quantities  across  the  wake.  The  mean  concentrations  and  the 
cylinder  siz<is  were  varied.  Injection  was  also  done  in  the  wake  of  the 
large  cylinder,  using  a  small  tube  as  a  "line  source". 

For  this  irreversible  reaction,  the  molecular  reaction  rate  term 
was  k  ,  where  k  was  a  constant  for  the  isothermal  experiment  and 

Cj  and  are  the  two  reactant  concentrations.  Thus  the  experimental 
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Table  in.  Chemical  Reactions  and  Rate 


3.5x  1012e-2550°/T 


situation  could  be  described  by  allowing  mean  and  fluctuating  concentra¬ 
tions  tc  be  averaged.  The  resulting  effective  reaction  rate  term  was 
k^l^2  +  <V<V).  w^ere  ^  and  ^2  are  mean  concentrations,  C^'  and 
C^'  are  the  fluctuating  concentrations  and  C^'C^7"  is  the  average  product 
or  correlation  coefficient.  An  equation  was  generated  for  1 1 ,  the 
terms  being  identified  and  assessed  as  to  their  significance.  An  effect 
of  incomplete  mixing  was  identified  that  later  proved  to  be  especially 
significant  for  calculations  of  fast  reactions.  For  incomplete  mixing, 

tended  to  become  negative,  thus  reducing  the  effective  reaction  rate 
k^l^2  +  below  which  is  the  rate  for  the  ideal  case  of  com¬ 

plete  mixing  on  a  molecular  level.  Unmixed  reactants  passing  by  a  point 
in  the  flow  field  contribute  to  the  mean  concentration  term,  but  the  actual 
reaction  was  only  taking  place  on  the  edge  of  the  "clumps"  of  the  reactants 
where  just  a  few  of  the  reactant  molecules  were  mixed.  It  should  be  noted 
that  the  effective  turbulent  reaction  rate  is  still  higher  than  that  for  a 
laminar  wake  flow  because  the  reactants  are  turbulently  dispersed  over 
a  larger  volume  and  partially  mixed,  creating  many  more  sharp  gradients 
for  molecular  diffusion  than  would  exist  in  a  laminar  flow.  The  compari¬ 
son  of  the  equations  with  the  theory  was  limited  because  only  one  of  the 
two  reactants  was  measured.  It  was  evident,  however,  that  the  effect  of 
incomplete  mixing  compared  with  ideal  complete  mixing  was  similar  to 
the  effect  of  a  nonequilibrium  reaction  in  that  the  reaction  is  delayed. 
Further,  it  also  is  evident  that  the  1  term  which  was  typical  of  the 

Rubel-Zeiberg  terms,  was  a  variable  throughout  the  flow  field  and  not  a 
constant. 

The  fast-reaction  theory  of  Toor  (Reference  39)  was  investigated 
toward  the  end  of  the  contract  period.  It  treats  the  limiting  case  of 
reactions  that  are  fast  compared  with  the  rate  of  molecular  diffusion. 

Parts  of  the  Toor  theory  were  developed  in  the  study  of  turbulent  flames 
(Reference  40)  which  are  an  example  of  diffusion-limited  reactions.  For 
isotropic  turbulent  mixing  that  was  induced  by  a  grid  in  ?,  pipe,  the  theory 
has  been  verified  (Reference  41)  experimentally.  It  should  be  noted  that 
the  theory  requires  mean  and  fluctuating  values  of  one  reactant  undergoing 
diffusion  or  mixing  without  the  other  being  present.  Thus,  a  mixing 
experiment  is  an  input  and  the  mean  concentrations  during  the  reaction  is 
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the  output  of  the  theory.  The  actual  kinetics  of  the  reaction  are  elimin¬ 
ated  from  the  problem  by  a  transformation  of  variables  that  is  based  on 
equilibrium  kinetics  for  the  reaction  rate  terms  in  the  instantaneous  spe¬ 
cies  conservation  equations.  Thus  the  actual  reaction  rate  need  not  be 
known  if  it  is  fast  compared  to  diffusion. 

The  Toor  theory  was  compared  with  the  results  of  the  cylinder-wake 
experiment  previously  discussed.  It  should  be  noted  that  the  molecular 
diffusion  and  reaction  rates  were  not  measured  in  this  experiment.  More¬ 
over,  they  were  not  found  in  aliterature  search  ccnducted  during  this 
study.  The  fact  that  the  reaction  is  one  of  the  fastest  known  to  chemical 
analysts  led  to  the  belief  that  the  experiment  may  be  diffusion-  controlled 
and  thus  be  comparable  to  the  theory  of  Toor.  The  major  agreement 
between  the  Toor  theory  and  the  cylinder -wake  experiment  occurred  on 
the  wake  axis.  There  the  effective  reaction  rate  from  both  the  measure¬ 
ments  and  theory  was  higher  because  the  mixing  was  greater. 

It  is  clear  from  the  low-speed  water-tunnel  studies  undertaken  thus 
far  that  diffusion  or  mixing  must  be  studied  at  the  same  time  as  fluctua¬ 
tion  effects  on  nonequilibrium  chemistry.  A  practical  calculation  Scheme 
for  low-speed  isothermal  flows  must  specify  the  relation  between  dif¬ 
fusion,  reaction  and  convection  for  each  species  and  reaction  occurring 
at  each  point  in  the  turbulent  flow.  Some  additional  progress  can  be  made 
by  continued  studies  of  the  Toor  theory  for  fast  reaction  and  the  Rubel- 
Zeiberg  theory  for  slow  reactions,  as  described  above.  It  is  clear,  how¬ 
ever,  that  a  complete  description  of  the  reentry-vehicle  boundary-layer 
flow  and  chemistry  awaits  additional  experimental  data  on  chemical 
reactions  in  gases  undergoing  turbulent  flow.  There,  the  effects  of  tem¬ 
perature  and  density  fluctuations  on  the  effective  reaction  rate  will  be 
revealed. 

4,4.3  Wind-Tunnel  Kinetics  Experiment 

Under  a  TRW -sponsored  program,  a  low  speed  wind  tunnel  has  been 
designed  and  is  in  the  process  of  being  acquired  with  TRW  funds.  Experi¬ 
ments  using  dissociation  and  recombination  of  nitrogen  oxides  in  turbulent 
flow  at  elevated  temperatures  are  being  planned.  Both  mean  and  fluctu¬ 
ating  concentrations  of  NO^  will  be  measured  locally  with  a  light  probe. 
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The  attenuation  of  light  at  the  proper  wavelength  passing  over  a  short  gap 
in  the  flow  will  be  a  direct  measure  of  the  NC^*  the  other  species  being 
transparent  at  that  wavelength.  The  light  path  in  the  flow  will  be  between 
the  ends  of  a  fiber  optics  probe.  Thus,  the  averaging  length  of  the  light 
beam  will  be  a  minimum  consistent  with  the  sensitivity  of  the  instrumen¬ 
tation.  Conditions  in  the  tunnel  will  be  established  so  that  a  wide  range 
of  reaction  rates  in  turbulent  flow  can  be  examined  by  utilizing  a  range  of 
elevated  temperatures  and  flow  configurations.  Buoyancy  effects  will  be 
minimized  by  utilizing  vertical  flow  in  the  test  section.  Thus,  the  needed 
experimental  data  on  fluctuation  effects  on  chemical  reactions  in  turbulent 
gaseous  flow  will  be  obtained. 

4.5  CONCLUSIONS 

Considerable  progress  has  been  made  in  each  of  the  areas  of  this 
task.  The  approach  originally  outlined  still  appears  to  offer  the  most 
promise  with  respect  to  accomplishing  the  overall  objective  of  this  task. 

4.6  RECOMMENDATIONS 

TRW  recommends  that  the  following  steps  be  taken  to  complete  the 
analysis: 

Include  the  effects  of  viscous  dissipation  and  intermittency 
in  the  low  speed  boundary  layer  analysis. 

Derive  the  actual  scaling  functions  necessary  to  model  a 
high  speed  flow  to  its  low  speed  equivalent. 

Apply  the  transformation  to  the  low  speed  analytical  model. 

Complete  the  analysis  required  for  the  determination  of  the 
chemical  kinetics. 

Include  chemistry  in  the  analytical  boundary  layer  model. 

Verify,  to  whatever  extent  it  is  feasible,  this  model  by 
comparison  to  experiment. 
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SECTION  V 


IMPROVED  LAMINAR  WAKE 

(HYPERSONIC  WAKE  AND  BOUNDARY  LAYER  TECHNOLOGY) 

5.1  INTRODUCTION  AND  BACKGROUND 

The  objective  of  this  task  is  to  obtain  a  more  accurate  descrip¬ 
tion  of  the  extended  slender  body  far  wake  in  order  to  interpret  recent 
ICBM  flight  data  and  ballistic  range  data,  as  related  to  the  "break¬ 
through"  problem  (i.e.  ,  the  radar  return  from  slender  body  wakes  seems 
to  imply  that  the  turbulent  inner  viscous  core  is  shielded  by  an  overdense 
rotational  outer  wake  for  a  considerable  distance  downstream  of  the 
transition  zone,  and  thus  the  turbulent  front  of  the  inner  core  is  not  ob¬ 
served  until  it  "breaks  through"  this  ionized  portion  of  the  outer  wake). 

The  accomplishment  of  the  Improved  Laminar  Wake  task  required 
two  major  activities: 

(a)  programming  of  the  TRW  finite  difference  wake 
solution  for  50  radial  mesh  points, 

(b)  prediction  of  "breakthrough"  and  comparison  with 
test  data. 

The  first  of  the  above  activities  was  initially  completed  for  the  G.  E. 

635  computer.  However,  systems  difficulties  associated  with  this  com¬ 
puter  so  degraded  its  reliability  that  the  wake  program  was  finally 
transferred  to  the  C.  D.  C.  6600  machine.  A  sample  calculation  has  been 
made  to  illustrate  the  scheme  for  describing  the  complete  slender  body 
far  wake,  and  the  details  of  this  method  are  present  in  5.  3. 

The  TRW  wake  program  of  interest  uses  a  finite  difference  method 
to  solve  the  viscous  wake  equations  for  nonequilibrium  flow.  Initial 
conditions  and  outer  edge  conditions  are  arbitrary  input  data,  and  either 
laminar  or  turbulent  transport  properties  may  be  used.  The  conservation 
equations,  transformations,  and  numerical  methods  utilized  in  this 
program  are  presented  in  Reference  42. 

The  laminar  version  of  the  nonequilibrium  wake  program  considers 
seven  air  species  {N2,  N,  0^,  NO,  NO+,  e  )  plus  one  inert  contaminant. 
Local  equilibrium  of  vibrational  and  electronic  excitation  is  assumed  for 
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all  species.  Ambipolar  diffusion  is  adopted  for  the  NO  -e  pair,  and 
the  complete  multicomponent  equations  of  Reference  43  are  used  to 
calculate  the  transport  properties  of  the  gas  mixture  (including  thermal 
diffusion).  The  turbulent  version  of  the  program  uses  transport  prop¬ 
erties  based  on  the  eddy  diffusivity  concept  (Reference  44)  and  also 
includes  the  two  additional  species  O  and  O  ^ • 

5.  2  TECHNICAL  APPROACH 

One  of  the  more  interesting  features  of  the  far  wake  is  that  the 
radial  pressure  gradient,  ^p/3r»  is  small  across  both  the  inner  and  outer 
wake.  This  implies  that  the  classical  boundary  layer  approximations  to 
the  Na vie r -Stokes  equations  can  be  used  for  both  flow  regions.  Hence, 
the  major  effort  within  this  task  was  directed  toward  the  problem  of  ex¬ 
tending  an  existing  TRW  computer  program  which  currently  describes 
the  inner  viscous  core  of  the  nonequilibrium  laminar  wake  downstream 
of  the  critical  point,  radially  outward  into  that  cold  gas  portion  of  the 
wake  which  is  both  upstream  of  transition  and  within  the  envelope  of 
the  wake  and  body  shocks  (Figure  49). 

The  techniques  discussed  in  Section  VI  were  employed  to  describe 
the  outer  flow  at  the  beginning  of  the  laminar  far  wake  region;  initial 
conditions  for  the  inner  viscous  core  will  be  estimated  by  using  test 
data.  The  TRW  integral  method  was  utilized  for  the  turbulent  wake. 

As  noted  above,  the  radial  extension  of  the  laminar  wake  solution 
does  represent  an  improvement  over  the  traditional  method  of  allowing 
a  viscous  region  to  grow  into  a  rotational  inviscid  flow  in  that  once  the 
initial  conditions  are  established,  no  arbitrary  separation  of  viscous 
laminar  andinviscid  rotational  flow  regimes  is  required;  boundary  con¬ 
ditions  are  applied  in  the  cold  undissociated  flow  at  the  extreme  edge  of 
the  outer  wake.  Thus,  the  problem  of  matching  gradients  at  the  inter¬ 
face  between  the  inner  and  outer  wake  is  avoided.  The  rate  ai  which 
the  inner  and  outer  wakes  merge  is  governed  by  the  molecular  trans¬ 
port  properties  of  the  gas  and  was  determined  in  the  course  of  the 
calculation. 
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Figure  49.  Flow  1-  \eld  Sketch 


5.3  RESULTS 


A  far  wake  calculation  was  performed  to  compare  "breakthrough" 

with  the  observed  ballistic  range  results  of  Reference  58.  The  body  was 

a  12.  5  half- angle  cone  of  .  06  bluntness  ratio  with  a  base  diameter  of 

.  25  inches.  A  velocity  cf  23,  400  ft/sec  and  an  ambient  pressure  of  75 

torr  were  considered.  The  freestream  Reynolds  number  R  was 
5  e, 

3  x  10  .  Transition  of  the  viscous  inner  core  from  laminar  to  turbulent 
was  assumed  to  occur  at  50  X/D. 

General  results  of  the  calculation  are  not  presented  herein;  how¬ 
ever,  the  radial  enthalpy  variation  at  the  beginning  and  end  of  the  lam¬ 
inar  wake  are  indicated  in  Figure  50. 

Growth  of  the  turbulent  core  is  shown  in  Figure  51  along  with  the 
predicted  inviscid  edge  boundary.  This  calculation  indicates  "break¬ 
through"  to  occur  about  300  body  diameters  downstream  of  the  base. 

This  is  in  agreement  with  the  distance  observed  in  Reference  58. 

5.4  CONCLUSION 

In  conclusion,  the  TRW  finite  difference  wake  program  has  been 
modified  to  provide  better  radial  definition  of  the  laminar  far  wake.  This 
program  was  successfully  utilized  in  a  calculation  to  determine  the  loca¬ 
tion  of  "breakthrough";  prediction  and  observation  agreed  satisfactorily. 
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Figure  50.  Radial  Enthalpy  Variation 
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SECTION  VI 


NEAR  WAKE 

(HYPERSONIC  WAKE  AND  BOUNDARY  LAYEF  TECHNOLOGY) 

6.  1  INTRODUCTION  AND  BACKGROUND 

The  near  wake  of  a  slender  body  in  hypersonic  flow  is  a  region  of 
strong  viscous  interaction  where  large  changes  in  pressure,  flow  direc¬ 
tion,  velocity,  and  temperature  occur  over  but  a  few  base  diameters.  It 
is  the  region  which  literally  "sets  the  stage"  for  the  flow  hundreds  of  base 
diameters  downstream,  where  much  of  the  interest  is  on  radar  observ¬ 
ables.  In  fact,  electrons  and  contaminants  from  the  boundary  layer  pass 
into  the  near  wake  where  their  concentrations  are  altered.  In  addition, 
electrons  themselves  can  be  produced  due  to  the  presence  of  local  high 
temperatures  brought  about  by  viscous  dissipation. 

Since  the  near  wake  is  instrumental  in  controlling  the  downstream 
flow,  it  must  be  a  part  of  the  development  of  a  mathematical  model  for 
wake  flows.  However,  the  presence  of  many  coupled  and  complex  fluid 
dynamic  phenomena  in  this  relatively  small  region  has  made  its  analysis 
very  difficult.  As  a  result,  various  aspects  have  been  studied  in  turn  at 
TRW  with  the  intent  of  establishing:  1)  the  basic  mathematical  methods 
which  can  be  used  in  the  overall  problem,  2)  models  of  various  portions 
of  the  flow  field,  and  3)  coupling  schemes  for  joining  the  various  regions. 
Background  and  more  detailed  descriptions  of  these  efforts  are  discussed 
in  the  following  paragraphs. 

In  the  near  wake  region  of  a  hypersonic  slender  body,  shown  in 
Figure  52  ,  we  are  concerned  with  the  analytical  determination  of  the 
steady  state  flow  field.  There  are  two  basic  approaches  to  this  problem. 
The  first  is  to  treat  the  steady  state  as  the  asymptotic  limit  of  a  time 
dependent  solution.  If  we  limit  ourselves  to  continuum  flow,  the  problem 
becomes  mathematically  parabolic  in  time  with  two  space  variables. 

There  are  some  conceptual  difficulties  with  the  determination  of  down¬ 
stream  boundary  conditions  but  this  approach  seems  attractive  for  prob¬ 
lems  with  a  very  small  Reynolds  number.  However,  for  the  Reynolds 
numbers  we  are  interested  in,  the  computing  speed  and  storage  capacity 
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needed  to  obtain  the  required  spatial  resolution  are  beyond  the  capabilities 
of  current  computers. 

The  second  approach  to  this  problem  is  to  formulate  it  directly  in 
terms  of  steady  state  equations.  We  will  concern  ourselves  here  only 
with  continuum  laminar  flow  described  adequately  by  the  steady  state 
Navier-Stokes  equations.  The  formal  mathematical  nature  of  these 
equations  indicates  that  this  then  becomes  a  boundary  value  problem;  that 
is,  its  solution  requires  the  specification  of  boundary  conditions  along 
all  boundaries.  If  no  further  simplifying  assumptions  are  made,  the 
problem  formulated  in  this  way  suffers  from  the  same  difficulties  at  the 
Reynolds  numbers  of  interest  as  does  the  time  dependent  approach.  For¬ 
tunately  however  there  are  regions  of  the  flow  in  which  some  of  the  higher 
order  terms  of  the  Navier-Stokes  equations  become  negligible.  For 
instance,  the  flow  over  the  major  portion  of  the  body  can  be  divided  into 
a  boundary  layer  and  an  outer  inviscid  supersonic  flow,  with  the  equations 
describing  each  region  obtained  by  neglecting  higher  order  terms  in  the 
Navier-Stokes  equations.  The  boundary  layer  equations  have  a  parabolic 
character,  the  inviscid  supersonic  flow  equations  are  hyperbolic,  and 
both  are,  of  course,  initial  value  problems.  If  the  requisite  boundary 
conditions  are  known  a  priori,  as  they  are  for  these  two  problems,  the 
calculation  proceeds  in  a  straightforward  manner  by  marching  in  a  gen¬ 
erally  downstream  direction. 

6.2  CORNER  FLOW 

As  the  body  shoulder  is  approached,  this  rather  simple  flow  model 
begins  to  become  inadequate  since  it  is  then  no  longer  possible  to  divide 
the  flow  into  a  constant  pressure  viscous  (boundary  layer)  region  and  an 
exterior  inviscid  flow.  The  upstream  influence  of  the  low  base  pressure 
on  the  approaching  boundary  layer  produces  a  strong  streamwise  pressure 
gradient  and  a  corresponding  cross-stream  pressure  gradient  in  the 
supersonic  viscous  region.  Considerable  progress  has  been  made  under 
the  current  Penetration  Systems  Studies  contract  toward  the  understanding 
of  this  flow  region  and  the  development  of  methods  of  computing  the  flow 
properties  in  detail.  Reference  45  describes  the  analysis  of  flows 
approaching  a  sharp  corner;  Reference  46  includes,  in  addition,  an 
analysis  of  the  flow  up  to  and  around  a  rounded  corner  to  the  separation 
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point.  In  these  analyses,  the  boundary  layer  equations  are  used  to 
describe  the  subsonic  portion  of  the  flow  with  the  effect  of  transverse  pres¬ 
sure  gradients  in  the  supersonic  viscous  region  included  through  the  use 
of  an  inviscid  transverse  momentum  equation.  A  relation  between  pres¬ 
sure  and  flow  angle  at  the  outer  edge  is  obtained  by  coupling  with  an  outer 
inviscid  region.  Initial  profiles  are  assumed  to  be  similar  (Blasius)  in 
shape,  but  with  the  boundary  layer  edge  flow  inclination  not  specified. 
When  the  equations  are  solved  by  an  implicit  finite  difference  method,  it 
is  found  that  the  qualitative  nature  of  the  downstream  flow  is  determined 
by  the  particular  value  of  initial  edge  flow  inclination  chosen,  with  two 
families  of  solutions:  either  there  is  an  unfavorable  pressure  gradient 
with  boundary  layer  separation  downstream,  or  there  is  a  precipitous 
drop  in  pressure  along  the  wall.  Bounding  these  two  families  of  solutions 
is  the  particular  solution  which  continues  downstream  on  the  body  without 
separation.  In  this  formulation,  higher  order  Navier-Stokes  terms  whose 
function  is  to  determine  the  location  of  the  separation  point  have  been 
neglected.  In  the  absence  of  these  terms,  solutions  can  be  generated 
corresponding  to  all  possible  locations  of  the  separation  point  and  the 
appropriate  solution  is  chosen  by  using  independent  information.  This  is 
illustrated  in  Figure  53  ,  which  shows  the  wall  pressure  distribution 
along  the  rounded  base  of  a  wedge.  The  experimental  pressures  meas¬ 
ured  by  Hama  (Reference  47)  are  given  by  the  bars,  and  the  solid  lines 
are  members  of  the  computed  families  of  pressure  distributions.  The 
family  for  which  the  pressure  reaches  a  minimum  and  then  begins  to  rise 
corresponds  to  solutions  separating  from  the  curved  wall  under  an 
adverse  pressure  gradient,  and  all  are  valid  solutions  corresponding  to 
different  values  of  base  pressure.  These  can  be  obtained  experimentally 
by  blowing  or  suction  .n  the  recirculation  region  to  change  the  base  pres¬ 
sure.  The  family  for  which  the  pressure  drops  rapidly  corresponds  to 
cases  in  which  the  wall  curvature  is  increased  just  downstream  of  where 
the  solution  is  terminated.  The  pressure  drop  corresponds  to  the 
upstream  influence  of  this  change  in  wall  curvature,  and  can  be  obtained 
experimentally  by  truncating  the  rounded  afterbody  and  rounding  the 
newly  created  sharp  corner. 

In  Figure  53  ,  the  solution  best  fitting  the  experimental  pressure 
distribution  after  the  separation  process  has  been  initiated  has  been 
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carried  out  to  the  separation  point  itself.  The  detailed  flow  field  corre¬ 
sponding  to  this  calculation  is  shown  in  Figure  54  ,  and  the  point  which 
is  immediately  apparent  is  that  the  calculation  produces  a  separation  (or 
lip)  shock  deeply  imbedded  in  the  boundary  layer.  In  Figure  53  ,  the 
computed  separation  shock  location  and  streamlines  are  superimposed 
over  a  shadowgraph  corresponding  to  the  same  flow  conditions.  The 
agreement  is  striking. 

A  further  examination  of  Figure  54  shows  that  viscous  effects  are 
important  out  to  about  that  streamline  which  is  initially  at  Mach  1.  This 
immediately  suggests  a  simplified  model  in  which  the  supersonic  bound¬ 
ary  layer  is  expanded  inviscidly  around  the  body,  say  by  using  the  method 
of  characteristics.  Unfortunately,  however,  an  examination  of  the  shape 
of  the  streamline  initially  at  Mach  1  discloses  that  there  is  no  obvious 
relation  between  its  shape  and  that  of  the  wall.  A  calculation  assuming 
that  this  streamline  follows  the  shape  of  the  wall  results  in  pressures 
which  are  far  below  those  measured.  It  therefore  appears  that  such  a 
simplified  model  is  inadequate. 

In  Figure  55  ,  the  computed  pressure  distribution  on  a  sharp 
cornered  wedge  is  compared  with  experimental  values  for  three  flow 
conditions,  again  with  strikingly  good  agreement, 

6.  3  DOWNSTREAM  EFFECTS  OF  UPSTREAM  INFLUENCE 

There  is  another  important  aspect  of  the  corner  flow  problem  which 
was  discovered  as  the  result  of  calculations  using  an  overly  simplified 
model  of  the  corner  expansion.  The  model  was  proposed  by  Weinbaum 
(®-e^erence  48),  and  assumes  that  the  supersonic  boundary  layer  under¬ 
goes  an  inviscid  centered  expansion  at  the  body  corner  from  the  body 
surface  pressure  to  the  base  pressure.  Such  a  model  ignores  the 
upstream  influence  cf  the  base  pressure  on  the  boundary  layer.  Calcu¬ 
lations  of  this  type  were  performed  at  TRW  using  the  method  of  rotational- 
characteristics  (Reference  49)  with  initial  conditions  which  were  consistent 
with  experiments  conducted  by  Batt  (Reference  50)  on  a  20-degree  wedge 
at  a  free  stream  Mach  number  of  6.  Batt  was  able  to  resolve  all  of  the 
details  of  the  flow  (static  and  total  pressure  and  temperature,  streamline 
orientation,  and  velocity  distribution)  in  the  near  wake  beginning  just 
downstream  of  the  corner.  Some  of  his  results  are  shown  in  Figure  56. 
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Figure  55.  Computed  Streamlines  and  Separation 
Shock  Superimposed  on  Shadowgraph 
Corresponding  to  the  Same  Conditions 
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Figure  56.  Wall  Pressure  Near  a  Sharp  Corner-Comparison 
of  Theory  with  Experiment 
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Comparisons  of  the  characteristics  calculation  with  the  experiment  over 
a  distance  of  several  base  heights  downstream  of  the  corner  were  made, 
and  the  results  indicated  that  the  centered  expansion  grossly  overpredicts 
(by  a  flow  angle  of  as  much  as  10  degrees)  the  streamline  orientation  in 
the  inviscid  region  of  the  flow  where  the  expansion  waves  from  the  corner 
have  an  influence.  In  Figu>  e  57  ,  this  inviscid  region  lies  above  the 
dotted  line  at  approximately  1/2  of  the  total  base  height  (Y/H  ~  0.  5)  above 
the  X/H  axis.  This  line  represents  the  locus  of  points  where  the  change 
in  total  pressure  along  streamlines,  hence  viscous  dissipation,  becomes 
important.  In  order  to  obtain  gcod  agreement  with  the  experimental  data, 
it  was  necessary  to  disregard  the  model  employing  the  centered  expan¬ 
sion,  and  in  its  place  consider  a  non-centered  expansion  beginning 
upstream  of  the  corner.  In  the  latter  calculation,  the  streamline  geom¬ 
etry  forming  the  inner  boundary  of  the  inviscid  flow  to  be  calculated  is 
specified  as  a  smooth  curve,  rather  than  a  straight  line  with  a  discontin¬ 
uous  slope  near  the  corner  (as  in  the  centered- expansion  case).  The  non- 
centered  expansion  reflects  the  upstream  influence  of  the  base  pressure 
on  the  boundary  layer  structure  before  the  corner,  and  also  demonstrates 
how  this,  in  turn,  has  an  important  effect  on  the  downstream  flow  enter¬ 
ing  the  viscous  near  wake. 

6.4  INVISCID  FLOW  ASSUMPTION 

Downstream  of  the  separation  point  on  the  body,  an  analytical  model 
of  the  near  wake  was  formulated  at  TRW  under  the  assumption  that  viscous 
forces  are  important  in  the  recirculation  region,  in  a  narrow  shear  layer 
adjacent  to  the  dividing  streamline,  and  below  the  wake  shock  wave.  The 
flow  above  these  regions  (see  Figure  52  )  was  assumed  to  be  expanding 
inviscidly;  the  rationale  for  this  is  that  the  vorticity,  which  was  initially 
very  large  in  the  boundary  layer,  would  be  significantly  and  rapidly 
decreased  during  the  expansion  to  the  point  where  viscous  dissipation 
would  be  unimportant.  As  a  result,  a  method-of-characteristics  program 
was  written  (Reference  49)  for  rotational  flow  which  included  the  capability 
of  calculating  imbedded  shock  waves.  An  example  showing  the  right 
running  characteristics  describing  the  expansion  of  the  boundary  layer 
into  the  near  wake  is  presented  in  Figure  58.  A  lip  shock  is  shown, 


igure  57.  Experimentally  Determined  Flow  F'eld 
in  the  Near  Wake  of  a  20-degree  Wedge; 
Mco  =  6,  Re  „  =  5,  5  x  10^  (Reference  5) 


and  a  wake  shock  results  as  a  consequence  of  an  imposed  pressure  rise 
along  the  inner  boundary  of  the  calculated  region. 

With  the  very  recent  experimental  data  of  Batt,  it  has  been  possible 
to  assess  the  validity  of  the  assumption  that  the  flow  is  expanding  invis- 
cidly.  It  was  stated  earlier  that  rotational  characteristics  calculations 
were  performed  with  initial  conditions  which  were  consistent  with  one  of 
Batt's  cases.  By  adjusting  the  curvature  of  the  bounding  streamline 
upstream  of  the  body  corner  all  flow  properties  at  a  level  of  a  base  height 
agreed  with  the  experiment.  The  shape  of  the  inner  edge  streamline  after 
the  corner  was  made  to  agree  with  the  experiment;  and  the  flow  field 
below  a  base  height,  and  above  where  the  wake  shock  would  be,  was  cal¬ 
culated.  The  result  is  presented  in  Figure  59,  where  the  right  running 
characteristics  are  given  and  two  isobars  are  shown;  the  upper  one  corre¬ 
sponds  to  the  base  pressure  and  the  lower  one  to  half  the  base  pressure. 
The  pressure  curve  shown  below  the  characteristics  represents  the 
calculated  pressure  along  the  specified  streamline  boundary  and,  when 
this  ends,  along  the  non- streamline  boundary  which  is  the  left  running 
characteristic  lying  just  above  the  wake  shock  (not  shown).  It  is  seen  that 
the  pressure  has  dropped  to  a  third  of  the  base  pressure  at  the  end  of  the 
calculation.  Thus,  the  flow  in  the  region  of  interest  is  expanding  and  this, 
unfortunately,  is  contrary  to  the  experimental  data.  The  Batt  data 
(Reference  50)  shows  that  the  pressure  is  constant  and  equal  to  the  base 
pressure  in  this  region,  and  that  there  are  significant  total  pressure 
losses  on  streamlines.  It  appears,  short  of  questioning  the  data,  that 
the  disagreement  can  only  be  explained  by  the  viscous  effects  which  are 
evident  from  the  data. 

The  experimental  data  suggests  that  the  viscous  near  wake  below  a 
base  radius  divides  itself  into  two  regions  in  which  pressure  forces  take 
on  different  relative  importance.  Above  the  wake  shock  the  balance  is 
between  shear  and  inertia  forces  whereas  the  region  below  the  shock  is 
where  the  primary  viscous  interaction,  and  hence,  pressure  rise  occurs. 
Evidently,  the  two  regions  are  mutually  coupled.  Clearly,  conditions 
above  the  shock  affect  those  below  it  since  the  flow  progresses  in  that 
direction.  On  the  other  hand,  the  presence  of  strong  viscous  dissipation 
in  the  upper  region  suggests  that  shearing  stresses  are  being  transmitted 
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Note:  Pressures  are  normalized  with  the  base  pressure 


across  the  shock  (probably,  in  a  direction  tangential  to  it).  Thus,  the 
basic  fluid  physics  of  the  near  wake  differs  from  that  assumed  with  the 
inviscid  flow  nodel  above  the  wake  shock.  In  that  model,  the  balance  is 
between  pressure  and  inertia  forces,  and  the  flow  above  the  shock  is 
influenced  by  left  running  characteristics  which  originate  from  the  edge 
of  the  shear  layer  above  the  dividin ;  streamline  within  close  proximity  to 
the  base,  and  right  running  characteristics  which  emanate  from  the  flow 
above  the  corner.  However,  the  experimental  results  show  that  this  is 
fundamentally  in  error.  Pressure  forces  are  not  important  above  the 
shock,  and  the  flow  there  is  strongly  influenced  by  the  flow  below  the 
shock.  It  will  be  shown  later  that  the  solution  to  the  interacting  near 
wake  is  very  sensitive  to  the  flow  conditions  above  the  shock,  and  hence, 
that  accurate  determination  of  the  flow  in  this  region  is  very  important  to 
the  overall  wake  problem.  It  must  be  concluded  that  attempts  at  a  near 
wake  solution  which  employ  the  inviscid  model,  such  as  the  partial  solu¬ 
tions  of  Weiss  (Reference  51),  are  therefore  of  questionable  validity. 

6.  5  WAKE  NECK 

The  flow  downstream  of  the  wake  stagnation  point  has  been  studied 
in  some  detail  at  TRW  and  is  relatively  well  understood.  A  flow  model 
coupling  a  boundary-layer-like  region  near  the  axis  to  an  inviscid  outer 
region  using  a  displacement  interaction  was  first  investigated  by  Webb, 
Golik,  Lees,  and  Vogenitz  (References  52,  53).  They  found  that  for  any 
one  parameter  family  of  initial  profiles,  the  solution  downstream  could 
be  divided  into  two  families,  with  a  particular  solution  forming  the  bound¬ 
ary  of  these  two  families  and  proceeding  downstream  through  a  saddle 
point  type  of  singularity.  Only  the  particular  solution  was  found  to  have 
the  downstream  behavior  appropriate  to  a  far  wake  and,  since  the  param¬ 
eter  describing  the  initial  profile  shapes  could  be  related  to  the  base 
pressure,  this  introduces  a  required  uniqueness  condition  which  deter¬ 
mines  the  base  pressure. 

This  behavior  is  not  at  all  surprising,  since  there  are  other  well- 
known  examples  of  problems  having  qualitatively  similar  solutions.  For 
instance,  the  one  dimensional  nozzle  problem  is  hyperbolic  when  the 
time  dependent  term  is  retained  and  the  steady  state  solution  is  the 
asymptotic  limit  for  large  time.  Whether  one  finds  a  supersonic  or 
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subsonic  solution  depends  on  the  imposed  boundary  condition  at  the  nozzle 
exit.  Howe  /er,  when  one  searches  for  a  solution  to  the  steady  state 
equations,  it  is  found  that  this  becomes  a  problem  with  parabolic  charac¬ 
ter,  with  all  possible  downstream  conditions  obtainable  by  varying  an 
otherwise  undetermined  initial  condition.  There  is  a  whole  family  of 
subsonic  solutions,  but  the  supersonic  solution  in  unique  and  determined 
by  a  saddle -point  type  singularity.  If  we  start  with  the  knowledge  that 
the  downstream  flow  is  supersonic,  then,  in  either  the  wake  or  nozzle 
problem  the  only  solution  of  interest  io  the  singular  one. 

The  mathematical  nature  of  the  singularity  in  the  wake  problem  as 
formulated  by  Webb  et  al.  has  been  studied  in  great  detail  by  Ai  (Refer¬ 
ence  54)  at  TRW. 

Baum  and  Denison  (Reference  55)  extended  the  interacting  wake 
model  to  include  transverse  pressure  gradients  in  the  supersonic  region 
with,  of  course,  no  qualitative  change  in  the  results.  This  model  (in 
which  the  equations  are  solved  by  finite  differences)  and  the  Webb  et  al. 
model  (which  uses  integral  methods)  have  both  been  coupled  (by 
Ohrenberger,  Reference  56)  at  the  outer  edge  to  a  rotational  character¬ 
istics  calculation  (Reference  49)  which  includes  provisions  for  internally 
generated  shocks.  It  was  originally  expected  that  the  rotational  charac¬ 
teristics  -  viscous  wake  model  would  suffice  as  an  acceptable  model  of  the 
near  wake.  However,  as  previously  discussed,  recent  experimental  data 
suggest  that  viscous  effects  are  not  negligible  in  a  region  lying  above  the 
wake  shock  and  within  the  region  for  which  we  have  used  an  inviscid  flow 
model.  Nevertheless  it  is  possible  to  use  this  model  in  a  restricted  way, 
and  in  conjunction  with  known  experimental  results,  to  investigate  the 
sensitivity  of  the  wake  neck  solution  to  the  flow  above  the  shock.  In 
addition,  comparisons  can  be  made  with  experiment. 

Rather  than  employing  rotational  characteristics  to  calculate  the 
entire  expansion  of  the  boundary  layer  into  the  near  wake,  the  approach 
is  to  apply  the  method  to  a  narrow  region  adjacent:  to  both  sides  of  the 
wake  shockwave.  The  fact  that  viscous  dissipation  is  present  above  the 
wake  shock  is  then  reflected  in  the  boundary  conditions  imposed  on  the 
upstream  side  of  the  inviscid  region.  The  assumption  is  that  although 
the  integrated  viscous  dissipation  is  important  for  the  entire  expansion 
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process,  the  amount  during  a  short  traverse  of  the  streamlines  is  small. 
Batt's  experimental  results  bear  this  out  since  the  divergence  between 
streamlines  (long  dashed  lines  in  Figure  57  )  and  lines  of  constant  total 
pressure  (short  dashed  lines  in  Figure  57  )  from  the  boundary  layer  to 
wake  shock  wave  is  large,  whereas  the  local  divergence  is  relatively 
small. 

The  approach  to  the  present  calculations  is  more  easily  explained  by 
referring  to  Figure  60.  From  the  experimental  data,  the  necessary  flow 
parameters  are  specified  along  a  left  running  characteristic  which  lies 
just  above  the  wake  shock  wave.  The  viscous  model  is  restricted  to 
calculating  the  flow  downstream  of  the  wake  stagnation  point.  However, 
part  of  the  interaction,  and  hence  pressure  rise,  occurs  forward  of  this 
point.  This  is  accounted  for  by  specifying  from  the  data,  the  pressure 
distribution  upstream  of  the  wake  stagnation  point  on  that  streamline 
which  enters  what  the  model  considers  to  be  the  edge  of  the  viscous  wake 
at  the  station  of  the  wake  stagnation  point.  Only  one  other  condition 
remains  and  that  is  specification  of  the  velocity  profile  at  the  wake  stag¬ 
nation  point.  The  chosen  profile  must  reflect  the  correct  mass  balance 
consistent  with  the  data.  The  downstream  interaction  is  then  calculated 
as  an  eigenvalue  problem  with  the  value  of  Reynolds  number  as  the 
unknown  quantity.  Results  are  shown  in  Figures  61  an'’  62  where  the  two 
families  of  non-wa.ce -like  solutions  discussed  previously  are  shown.  The 
Reynolds  number  of  the  case  in  Figure  61  is  smaller  than  that  in  Figure 
62  ;  the  lower  Reynolds  number  case  (Re^  =  50,  000)  shows  a  source 
type  solution  (the  centerline  velocity  eventually  decreases  downstream) 
which  would  yield  a  second  stagnation  point  downstream  of  the  wake  stag¬ 
nation  point.  Figure  61  »  representing  the  same  flow  at  a  higher  Reynolds 
number  (Re^  =  55,  000),  shows  the  sink  type  solution  in  which  the  pres¬ 
sure  begins  to  decrease  along  the  axis.  The  solution  Reynolds  number  is 
therefore  bracketed  by  these  twn  cases.  The  Reynolds  number  corre¬ 
sponding  to  the  experimental  case  is  55,  000  and  so  the  agreement  is 
quite  favorable.  It  is  also  of  interest  to  compare  the  calculated  pressure 
distribution  along  the  axis  with  the  data.  This  is  done  in  Figure  63  and 
shows  that  the  calculated  distribution  is  generally  on  the  high  side  of  the 
experimental  data  but  is  nonetheless-acceptable. 
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Figure  60.  Wake  Neck  Interaction  Analysis 
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iteracting  Wake  Neck  Calculation  Showing  a  Source 
ype  Solution  (Re  =  50,000) 


Figure  62.  Interacting  Wake  Neck  Calculation  Showing  a  Sink 
Type  Solution  (Re  ^  =  55,  000) 


AXIS  PRESSURE 
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Figure  63.  Results  of  Wake  Neck  Interaction  and  Comparison 
with  Experiment 
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Another  point  of  interest  1-  the  sensitivity  of  the  solution  to  the 
input  flow  parameters.  There  is  some  flexibility  in  the  specification  of 
the  input  data  due  to  the  uncertainties  in  the  experimental  data.  This  is 
particularly  true  along  the  left  running  characteristic  above  the  wake 
shock.  Three  cases  were  considered  and  the  input  and  results  are  sum¬ 
marized  in  Figure  64.  The  best  agreement  with  experiment  is  for  Case 

I  which  has  already  been  presented.  The  low  angle  distribution  for  Case 

II  is  the  same  as  for  Case  I  but  the  Mach  number  although  initially  the 
same  was  increased  to  a  maximu'1*  of  about  8  percent.  The  input  left 
running  characteristic  is  thereby  inclined  at  a  somewhat  smaller  angle 
to  the  horizontal.  The  perturbation  in  Mach  number  led  to  a  solution 
Reynolds  which  lies  between  43  and  50  x  10  which  is  slightly  lower  than 
that  found  in  Cas  2  I.  The  axis  pressure  is  also  higher  and  has  diverged 
significantly  from  the  experimental  values  (the  solution  pressure  distri¬ 
bution  is  bracketed  by  the  curves  shown).  For  Case  III  the  flow  angle 
was  decreased  2.  5  degrees  and  the  Mach  number  decreased  by  a  maxi¬ 
mum  of  4  percent.  The  left  running  characteristic  is  approximately  ^ 
horizontal.  The  solution  Reynolds  number  lies  between  32  and  38  X  10 
which  is  in  considerable  disagreement  with  the  experimental  value  of 

55,  000.  The  axis  pressure  is  seen  to  rise  much  too  steeply  compared  to 
experiment.  These  results  demonstrate  that  the  solution  is  very  sensitive 
to  the  flow  above  the  shock  wave.  Changes  in  Mach  number  of  a  few 
percent  and  in  flow  angle  of  2.  5  degrees  have  a  dramatic  effect  on  the 
solution.  As  previously  stated,  it  has  become  apparent  that  an  accurate 
description  of  the  flow  above  the  wake  shock  is  required  if  meaningful 
near  wake  solutions  are  to  be  obtained. 

6.6  RECIRCULATION  REGION 

Recent  theoretical  investigations  of  the  flow  in  the  recirculation 
region  have  been  concerned  with  finding  exact  numerical  or  approximate 
series  solutions  to  the  Navier-Stokes  equations.  The  major  objection  to 
these  analyses  lies  with  their  use  as  a  component  part  of  an  overall 
analysis  of  the  near  wake,  and  not  as  a  description  of  the  recirculation 
region  itself.  In  the  overall  near  wake  or  "coupled"  problem,  the  flow 
adjacent  to  the  recirculation  region  may  be  described  by  a  set  of  parabolic 
differential  equations  which  would  normally  be  solved  using  a  forward 
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Figure  64.  Sensitivity  of  Wake  Neck  Solution  to  the  Flow 
Above  the  Wake  Shock 


marching  scheme.  However,  obtaining  an  exact  numerical  solution  for 
the  recirculation  region  is  an  elliptic  problem  and  this  fact  precludes  the 
use  of  marcning  methods.  It  is  therefore  not  possible  to  simultaneously 
calculate  the  inner  recirculating  and  outer  shear  flows.  In  this  method 
the  nume  i_al  procedure  must  involve  a  guess  of  the  boundary  conditions 
for  the  elliptic  region  followed  by  a  calculation  of  the  external  parabolic 
flow.  Successive  guesses  or  iterations  are  required  until  the  unmatched 
flow  properties  along  the  entire  boundary  between  the  two  regions  agree 
to  within  some  desired  accuracy.  Whether  an  acceptable  convergence  of 
the  flow  properties  along  the  entire  dividing  streamline  can  be  achieved 
is  questionable.  Certainly,  such  a  scheme  is  very  inefficient  from  the 
standpoint  of  computer  time.  This  is  further  compounded  when  a  solution 
to  the  overall  near  wake  problem  is  sought  since  the  iteration  described 
above  is  subordinate  to  the  iteration  for  the  flow  which  proceeds  through 
the  critical  point. 

Recently,  a  simplified  analysis  of  the  recirculation  region  was 
presented  by  Reeves  and  Buss  (Reference  57)  in  which  the  streamfunction 
was  described  by  a  double  series  expansion  in  the  axial  and  lateral 
co-ordinate  about  the  wake  stagnation  point  truncated  by  terms  of  the 
fifth  power.  Relations  between  the  coefficients  of  the  series  were  found 
by  satisfying  the  boundary  conditions  and  the  Navier-Stokes  equations. 
However,  it  remained  that  several  of  the  coefficients  must  be  specified. 

A  study  of  the  use  of  this  method  for  the  recirculation  region  in  the  over- 

4  , 

all  near  wake  problem  revealed  that  it  suffers  from  the  same  matching 
difficulties  as  in  the  numerical  case.  That  is,  guesses  must  be  made  for 
the  unknown  coefficients  followed  by  a  calculation  of  the  outer  parabolic 
flow.  This  is  repeated  until  agreement  of  the  unmatched  flow  properties 
occurs  on  the  boundary  between  the  inner  and  outer  regions.  The  ques¬ 
tions  of  convergence  and  excessive  computer  time  again  arise. 

What  is  needed  is  a  simple  analysis  of  the  recirculation  region 
which  is  compatible  with  the  marching  scheme  for  the  outer  shear  flow 
and  which  does  not  require  sub-iterations  for  matching  the  flow  properties 
on  the  dividing  streamline.  Such  an  analysis  was  initiated  during  the 
present  study. 
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The  present  approach  employs  integral  methods  applied  to  the 
Navier -Stokes  equations.  The  nature  of  the  mathematical  problem  is 
thereby  changed  from  elliptic  to  initial  valued,  thus  permitting  the  use 
of  marching  methods  in  conjunction  with  the  external  flow.  In  addition, 
a  simultaneous  matching  of  the  conditions  (flow  direction,  velocity, 
shear,  and  pressure)  on  the  dividing  streamline  with  the  external  shear 
flow  may  be  accomplished.  In  this  approach  the  only  iterations  required 
in  the  overall  near  wake  problem  will  be  on  the  parameter  which  dictates 
whether  the  flow  proceeds  through  the  critical  point  downstream  of  the 
wake  stagnation  point. 

To  facilitate  use  of  the  integral  method,  several  important 
assumptions  are  made.  The  first  has  to  do  with  the  form  of  the  stream- 
function,  which  in  an  integral  method,  must  be  specified.  It  was  noted 
in  the  analysis  by  Reeves  and  Buss  that  a  streamline  pattern  character¬ 
istic  of  the  recirculation  region  can  be  described  by  a  polynomial  of 
relatively  low  order.  In  particular,  a  fifth  order  polynomial  of  odd 
powers  in  the  lateral  co-ordinate  was  sufficient.  Exploiting  this  result 
for  the  integral  analysis  leads  to  the  following  form  for  ijj  : 

+  (x.q)  =  4jj(x)  n+  vj^2(x)n3i|j+  ^3(x)n5 

where  q  is  a  dimensionless  lateral  co-ordinate  y/6  (x),  and  6(x)  is  the 
height  of  the  dividing  streamline  from  the  x-axis.  The  functions 
6(x),  ^j(x),  4j2(x),  and  ^3(x)  are  unknown  dependent  variables  in  the 
integral  method.  Since  ^(x,q)  is  zero  when  q  is  unity  for  all  x,  it 
follows  that 


^3(x)  =  -  +  X(x)  -  4>2(x) 

and  the  number  of  dependent  variables  are  thereby  reduced  by  one. 

Also,  ^(x)  and  can  be  represented  in  terms  of  6(x),  the  centerline 

velocity,  uq(x),  and  the  horizontal  component  of  velocity  on  the 


)j( 

In  the  approach  of  Reeves  and  Buss,  +1,  ^2»  and  ^3  are  truncated 
series  expansions  in  x  with  undetermined  coefficients,  and  q  is  the 
lateral  co-ordinate  y. 
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dividing  streamline,  u^(x),  and  can  take  on  more  physical  meaning  as 
follows: 

4/j(x)  =  6  (x)  uq  (x) 

and 


4>2(x)  =  -26  (x)  uq(x)  -  —  6(x)  ufi(x) 


6(x). 

The  governing  equations  for  the  three  unknown  variables  are 
obtained  by  integrating  the  x-momentum  equation  (Navier-Stokes),  and 
a  moment  of  this  equation  across  the  recirculation  region.  The  additional 
equation  is  given  by  matching  the  shear  on  the  dividing  streamline  with 
that  from  the  outer  flow.  However,  in  the  absence  of  the  actual  coupling, 
the  external  shear  must  be  a  specified  function  of  the  axial  co-ordinate. 
The  other  flow  quantities  on  the  dividing  streamline  (pressure,  velocity, 
and  direction)  are  matched  as  a  consequence  of  specifying  the  boundary 
conditions  on  the  external  shear  flow  in  the  coupled  case. 

The  final  assumption  to  be  imposed  is  that  the  pressure  is  inde¬ 
pendent  of  the  lateral  co-ordinate  across  the  recirculation  region,  i.  e.  , 
p  =  p(x).  From  the  Reeves  and  Buss  solutions  at  the  higher  Reynolds 
numbers,  this  appears  to  be  a  good  assumption  away  from  the  base.  A 
significant  lateral  pressure  variation  does  appear  on  and  near  the  base; 
however,  it  may  be  shown  that  the  average  value  of  9p/9x  (which  appears 
in  the  integral  equations)  across  the  recirculation  region  is  very  small 
and  is,  in  fact,  exactly  zero  on  the  base.  Consequently,  the  average 
axial  pressure  gradient  term  will  be  ignored,  initially;  and  the  lateral 
pressure  field  near  the  base  will  be  calculated  using  the  y-momentum 
equation  once  the  solution  for  ^  is  known. 

The  governing  equations  have  been  developed  and  the  nature  of  the 
near  base  solution  investigated  by  a  series  expansion.  At  this  point,  the 
method  remains  promising  and,  hopefully,  will  be  continued  at  a  later 
date. 
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6.7  CONCLUSIONS  AND  RECOMMENDATIONS 

A  number  of  phenomena  have  been  identified  in  the  current  con- 
tract  period  an  being  of  critical  importance  in  determining  the  near  wake 
flow  field  These  include  two  major  effects: 


1) 

2) 


the  upstream  influence  of  the  low  base  pressure  on 
the  boundary  layer  approaching  the  separation  point 
strongly  affects  the  flow  field  far  downstream; 

viscous  effects  cannot  be  neglected  external  to  the 
lip  and  wake  shocks. 


A  method  of  calculating  the  boundary  layer  flow  approaching 
separation  has  been  developed,  including  the  upstream  influence  and  the 
lip  or  separation  shock  imbedded  deeply  within  the  boundary  layer.  The 
agreement  with  experiment  of  the  computed  wall  pressure  distribution 
near  separation  and  the  separation  shock  location  was  excellent.  It  is 
felt  that  the  phenomena  important  in  determining  the  near  wake  flow 
field  are  now  recognized,  and  the  tools  for  incorporating  them  into  a 
theoretical  model  have  been  developed.  The  most  important  problem  to 
be  attacked  in  the  near  future  is,  then,  to  tie  together  the  near  wake 
regions  into  a  complete  theoretical  near  wake  model  incorporating  all  of 
the  effects  which  we  now  recognize  to  be  important.  With  this  in  mind, 
the  investigation  of  a  simplified  method  of  computing  the  recirculation 
region,  suitable  for  incorporation  into  a  complete  near  wake  computa¬ 
tion  scheme,  has  been  initiated  and  is  continuing. 
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APPENDIX  II 


MONTE  CARLO  COMPUTER  SIMU L AT TON  OF  TWO-DIMENSIONAL 
AND  AXI-SYMMETRIC  SUPERSONIC  STEADY  FLOWS 


The  approach  is  to  conduct  numerical  experiments  with  a  model  gas 
on  the  computer.  The  real  gas  is  simulated  by  the  order  of  a  thousand 
rigid-sphere  molecules  which  may  be  thought  of  as  a  representative  sam¬ 
ple  of  the  many  billions  of  molecules  in  the  corresponding  real  gas.  The 
positions  and  velocity  components  of  the  simulated  molecules  are  stored 
in  the  computer  and  typical  collisions  are  computed  among  them  as  a  time 
parameter  is  advanced.  The  computation  of  collisions  starts  at  zero 
time,  the  molecules  having  been  set  up  as  a  uniform  stream  at  the 
required  freestream  Mach  number.  The  body  is  inserted  into  this  flow  at 
zero  time  and  the  desired  steady  flow  is  obtained  as  the  large  time  solu¬ 
tion  of  the  resulting  unsteady  flow. 

The  freestream  flow  is  in  the  positive  x  direction  and  the  simulated 
region  is  bounded  by  the  x  axis  as  a  line  of  symmetry,  and  by  the 
upstream,  outer  and  downstream  boundaries .  These  boundaries  must  be 
set  sufficiently  far  from  the  body  to  eliminate  interference.  The  simu¬ 
lated  region  is  divided  into  a  number  of  cells  which  are  sufficiently  small 
for  the  change  in  flow  properties  across  the  cell  to  be  small.  Since  the 
flow  is  either  two-dimensional  or  axi- symmetric*  only  two  position  coor¬ 
dinates  need  be  stored  for  each  simulated  molecule.  The  three  velocity 
components  are  also  stored  and  a  record  is  kept  of  the  molecules  within 

each  cell. 

The  first  step  is  to  generate  the  initial,  or  zero  time,  configuration 
of  molecules.  The  molecules  are  distributed  uniformly  over  the  simu¬ 
lated  region  and  the  velocity  components  are  appropriate  to  a  gas  in  Max¬ 
wellian  equilibrium  and  moving  at  the  required  Mach  number. 

The  body  is  then  inserted  into  the  flow  and  the  molecules  are 
allowed  to  move  and  collide  among  themselves.  The  two  processes  are 
uncoupled  by  computing  collisions  appropriate  to  a  time  interval  Atm  and 
by  then  moving  the  molecules  through  distances  appropriate  to  Atm  anc* 
their  instantaneous  velocities.  The  distortion  produced  in  the  molecular 
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paths  by  this  approximation  is  small  as  long  as  Atm  is  small  compared 
with  the  mean  time  between  collisions. 

Since  the  change  in  flow  properties  over  the  width  of  each  cell  is 
small,  the  molecules  in  a  cell  at  any  instant  may  be  regarded  as  a  sample 
of  the  molecules  at  the  location  of  the  cell.  The  relative  location  of  the 
various  molecules  within  the  cell  can  then  be  disregarded  and  the  colli¬ 
sion  probability  of  a  particular  pair  of  molecules  within  the  cell  depends 
only  on  their  relative  velocity.  A  pair  of  molecules  is  chosen  at  random 
from  those  within  the  cell  under  consideration  and  is  retained  or  rejected 
in  such  a  way  that  the  probability  of  retention  is  proportional  to  the  rela¬ 
tive  velocity  vr.  When  a  pair  is  retained,  a  typical  collision  is  computed 
between  the  two  molecules  and  the  new  velocity  components  are  stored  in 
place  of  the  old  ones.  The  random  selection  of  impact  parameters  is 
particularly  simple  for  rigid  sphere  molecules  since  all  directions  for  tb 
new  relative  velocity  vector  are  equally  probable. 

The  time  between  collisions  for  the  pair  is  (A  N  vr)  ,  where  N  is 
the  number  density  in  the  cell.  For  each  collision,  the  time  counter  for 

this  cell  is  advanced  by 

At  =  (2/Nc)  x  (A  Nvr)_1 


where  Nc  is  the  actual  number  of  molecules  in  the  cell.  Collisions  are 
computed  in  the  cell  until  the  time  counter  has  advanced  through  Atm. 
When  this  procedure  has  been  carried  out  for  every  cell,  the  overall  time 
is  advanced  through ■  Atm  and  the  molecules  are  moved  through  appropri¬ 
ate  distances.  Nc/2  is  of  course  the  number  of  pairs  in  the  cell.  Over  a 
large  number  of  collisions,  say  k,  the  average  time  between  collisions  is: 


At 


*  j-  1  <At>i  *  k*  I 


1 


i=l 


N  N4  Vr 
1  Ci  1  ri 


Now: 


1 


1  — k 


c.  "Ivr.  N  Nv 

l  l  c  r 


il 


k  N  N  v 

C _ £_ 


N  N,  v 
1  Ci  1  ri 


The  quantity  in  brackets  should  —  1  for  k  large  enough. 
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Then 


At  =  2/n  ~~m~  =  Tzr — \ 

c  ANvr  (nc/2) 

where  t  =  mean  time  between  collisions  per  molecule. 

Then  in  a  time  interval  "t,  on  the  average,  each  molecule  will 
undergo  one  collision  and  will  travel  an  average  distance  between  colli¬ 
sions 

t  •  c  =  — zrzr  ~  — ■ — 3  =  \ 

A  N  v  A  N 

r 

which  is  the  local  mean  free  path.  Note  that  by  normalizing  all  quantities 
to  a  reference  cell  in  the  free  stream,  the  molecular  cross  section  A 
need  never  be  specified  e.g., 

-A  _  ils  At  Jl  *^2  N<°  c=° 

*  N  ‘  =  Nc  N 

The  freestream  mean  free  path  is  determined  by  specifying  the  Knudsen 
numbe  r 


n  body  dimension 

The  set  of  molecules  in  each  cell  will  change  as  the  molecules  are 
moved  and  appropriate  conditions  must  be  applied  at  the  boundaries  of  the 
region  being  simulate  The  upstream  boundary  normal  to  the  free¬ 
stream  direction  is  treated  as  a  source  of  molecules  with  velocity  compo¬ 
nents  representative  of  the  downstream  moving  molecules  in  the  equilib¬ 
rium  freestream.  Any  molecule  which  moves  back  upstream  across  this 
boundary  is  regarded  as  being  "lost"  and  is  removed  from  the  store.  The 
plane  of  symmetry  along  the  x  axis  is  regarded  as  a  specularly  reflecting 
surface  in  two-dimensional  flows.  The  outer  and  rear  surfaces  present 
greater  difficulties.  A  procedure  has  been  developed  which  usually  intro¬ 
duces  only  a  smau  reflected  disturbance  and  which  becomes  exact  in  free 
molecule  flow.  This  is  to  regard  a  molecule  as  "lost"  if  it  moves  out- 
wa*"d  across  the  boundary  but,  if  it  moves  inward  from  the  boundary  dur¬ 
ing  the  time  interval  At^  through  a  distance  greater  than  its  original 
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distance  from  the  boundary,  it  moves  to  the  new  position  and  a  similar 
molecule  is  added  in  the  original  position. 

The  rationale  for  this  procedure  is  as  follows:  since  there  are  sup¬ 
posed  to  be  no  gradients  in  mean  properties  normal  to  the  boundary  (uni¬ 
form  flow)  it  is  equally  likely  that  a  similar  molecule  crosses  any  hori¬ 
zontal  line  near  the  boundary.  A  molecule  moving  into  the  field  away 
from  the  boundary  through  a  distance  greater  than  its  initial  distance 
from  it  is  a  statistical  measure  of  the  event  that  a  molecule  crosses  the 
boundary  and  enters  the  field.  Placing  it  so  that  it  enters  during  the  next 
movement  interval  is  a  convenient  way  of  handling  the  book  keeping  of 
getting  it  inside  the  field. 

The  latter  part  of  this  procedure  should  only  be  applied  to  those 
molecules  with  properties  similar  to  the  freestream  molecules.  An 
alternative  procedure  would  be  to  regard  the  outer  boundary  as  a  specu¬ 
larly  reflecting  surface.  While  this  would  result  in  a  larger  reflected 
disturbance,  it  does  permit  a  simple  physical  interpretation. 

Interactions  with  the  body  must  also  be  computed.  The  body  is 
assumed  to  have  a  diffusely  reflecting  surface  with  an  accommodation 
coefficient  of  unity.  After  the  flow  has  settled  down  to  a  steady  state,  the 
momentum  and  energy  transfer  to  the  surface  is  recorded  and  is  used  to 
compute  the  aerodynamic  data.  The  time  required  to  establish  steady 
flow  may  usually  be  assumed  to  be  close  to  the  time  required  for  the 
freestream  to  travel  several  body  dimensions.  The  flow  field  properties 
are  also  sampled  in  the  steady  flow.  Instantaneous  samples  are  recorded 
at  appropriate  time  intervals  and  these  are  averaged  for  greater  accu¬ 
racy.  The  time  interval  for  sampling  Atg  should  be  of  the  order  of  the 
time  required  for  the  flow  to  traverse  one  cell. 

The  results  become  progressively  more  accurate  as  the  sample 
size  is  increased  and  the  statistical  scatter  is  reduced.  An  indication  of 
the  order  of  this  scatter  is  given  by  the  normal  distribution  result  that  the 
standard  deviation  is  equal  to  the  reciprocal  of  the  square  root  of  the 
number  in  the  sample,  68,3 %  of  the  results  should  lie  within  one  stand¬ 
ard  deviation,  95.5%  within  two  and  99.7%  within  three.  Successful  runs 
have  been  made  with  as  few  as  six  simulated  molecules  in  each  cell  and 
the  statistical  fluctuations  do  not  appear  to  induce  instabilities  in  the  flow. 
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The  key  to  the  simulation  is  the  procedure  which  selects  the  mole¬ 
cules  for  a  collision  pair  according  to  the  correct  probability  and  then 
advances  the  time  parameter  by  the  appropriate  amount  for  each  colli¬ 
sion.  In  this  wa/,  an  appropriate  set  of  collisions  is  computed  for  any 
distribution  function  of  molecular  velocities.  Apart  from  the  boundary 
conditions,  only  two  approximations  are  made.  The  first  is  the  uncou¬ 
pling  of  the  molecular  collisions  and  motion,  and  this  becomes  exact  as 
Atm  tends  to  zero.  The  t-  cond  is  the  assumption  of  uniform  properties 
over  the  cell  and  this  disappears  as  the  cell  size  tends  to  zero. 

Programs  have  been  developed  for  supersonic  flow  past, 

i)  flat  plate  at  zero  incidence, 

ii)  circular  cylinder  with  axis  normal  to  flow, 

iii)  sphere, 

iv)  sharp  nosed  wedge  at  zero  incidence 

and 

v)  sharp  nosed  cone  at  zero  incidence. 

The  overall  number  flux^dragj and  heat  transfer  coefficient  are 
determined  for  each  body,  together  with  their  distribution  along  or  around 
the  surface  of  the  body.  Flow  field  properties  are  sampled  and  number 
density,  flow  velocity  and  temperature  are  output  at  a  number  of  points  in 
the  flow.  As  an  in  Vication  of  the  degree  of  non-equilibrium,  the  tempera¬ 
tures  based  on  the  individual  velocity  components  are  also  output.  In 
addition,  the  velocity  distribution  may  be  sampled  at  a  number  of  points 
in  the  flow.  Simulated  equilibrium  temperature  probes  may  be  set  in  the 
flat  plate  and  cylinder  flows. 
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APPENDIX  III 


PARAMETERS  OF  THE  MONTE  CARLO  CALCULATION 

The  steady  state  flow  field  properties  of  a  body  will  depend  upon  a 
set  of  fluid  mechanical  parameters;  viz, 

Cj-j  =  f  (Kn>  M,  T^/Tro,  Surface  reflection  properties. 

Fluid  properties,  Body  shape,  etc.) 

The  Monte  Carlo  estimate  of  the  same  property,  denoted  by  f°r  given 

fixed  values  of  the  above  fluid  mechanical  properties  will  depend  upon  a 
set  of  calculation  parameters;  viz, 

Cp  =  f  (Sample  size,  Chance,  Movement  time.  Cell  size, 

Initial  Field  accuracy.  Field  Size) 

The  first  two  quantities  reflect  the  fact  that  the  calculation  is  a  statistical 
one.  The  sample  size  is  the  number  of  molecules  to  contribute  informa¬ 
tion  regarding  the  particular  property  to  be  determined,  such  as  density 
or  drag.  Computing  the  arithmetic  mean  of  the  sample  (the  information 
from  each  molecule  is  weighted  equally)  we  can  appeal  to  two  laws  of 
statistics1^  to  tell  us  something  about  how  the  sample  mean  should 
behave.  The  Strong  Law  of  Large  Numbers  says  that  the  sample  mean 
will  tend  toward  the  true  mean  as  the  sample  size  tends  to  infinity: 

Cn  - *►  Cn  as  SS  - -  « 

USS 

The  Central  Limit  Theorem  further  says  that  as  the  sample  size  tends  to 
infinity  the  distribution  of  the  sample  mean  will  tend  toward  a  Gaussian 
distribution  with  variance 

2 

2  °PP 


where  o'pp  is  the  variance  of  the  parent  population  from  which  the  sam¬ 
ples  are  taken.  Note  that  the  sample  mean  is  a  distributed  variable  which 
means  that  if  the  calculation  is  repeated  with  different  random  numbers 
but  fixed  values  of  all  parameters  including  sample  size,  in  general  the 
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sample  me.*fi  vvili  be  different.  The  mention  of  infinity  in  the  two  laws 
need  not  cause  concern;  it  is  surprising  how  small  a  sample  is  needed  for 
the  laws  to  take  hold. 

In  general  it  is  difficult  to  make  a  priori  estimates  of  the  dispersion 
of  the  sample  mean  since  the  dispersion  of  the  parent  population  is 
unknown. 

For  the  following  simple  situation  an  approximate  numerical  esti¬ 
mate  can  oe  made;  consider  ner  r-free  molecule  flow  past  a  flat  plate. 
Assume  diffuse  reflection  from  plate. 


X 


Molecules  in  the  flow  field  can  be  considered  nearly  in  Maxwellian 
equilibrium  with  distribution  of  U  velocity 


7^u-"u->2 


f(U)  = 


A  °u 


where 


CO 

"u  =  \[2 

Because  the  reflection  is  diffuse  the  tangential  momentum  transferred  to 
the  plate  by  each  molecule  which  strikes  it  is  mU. 


Average  number  flux  to  unit  area  of  plate 


cm  P  vm_ 

CO  CO  CO 

4m  ~  2  \(tt  m 

In  a  small  time  interval,  At,  then  the  average  number  of  molecules  to  hit 
unit  area  of  the  plate  is: 

PVfnco 

sample  size  SS  =  — ^ — r - At 

2  «s/tt  m 
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The  estimated  momentum  flux  to  the  plate  is  then 


SS 


Momentum  flux 


5  *  it  I  m  ui  * 


PooVmoc  I 
2  sJtt  SS 


SS 


2  ui 


i  =  l 


i=l 


Normalizing  to  get  an  estimate  of  drag  coefficient: 

SS 
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'D  ”  i  T.  2 
2  P  Ucc 


mco  l  V 
-=-?«§)  2  ui 
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Then: 


Variance 
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An  approximation  for  the  mean  value  of  gg  Z  U.  is  U^.  Substituting  in  the 
expression  for  we  get 

^  1 


'Dr  1  ,  's/tt  S 

free  molecule 


Then 


1 _ 1_ 

•^2S'/ii 


For  more  complicated  flows  where  the  molecules  in  the  field  have 
widely  differing  distributions,  no  such  numerical  estimate  a.  f.  made, 
but  <r  would  be  expected  to  be  larger.  On  purely  empirical  grounds  it  has 
be<  assumed  in  the  analysis  in  this  report  that: 


<r 
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D 


1 

/ss 
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where  SS  is  the  sample  size.  Judging  by  the  scatter  in  repeated  calcula¬ 
tions,  this  somewhat  overestimates  <r. 

We  now  return  to  the  other  parameters  which  influence  the  flow 
properties;  viz,  Movement  time,  Cell  size.  Initial  Field  Accuracy,  and 
Field  Size.  In  order  for  the  calculation  method  to  be  easily  utilized,  it 
would  be  hoped  that  the  results  not  be  extremely  sensitive  to  these 
parameter  ? .  This  indeed  turns  out  to  be  the  case.  Movement  time.  At 

m 

and  cell  size  are  somewhat  akin  to  the  stepsize  in  a  finite  difference  cal¬ 
culation  —  they  need  only  be  reduced  in  size  until  the  change  in  the  renult 
falls  below  a  certain  level  (unlike  the  effect  of  step  size  the  calculation 
does  not  diverge  suddenly  as  they  are  increased,  the  degradation  of  accu¬ 
racy  occurs  rather  slowly). 

The  movement  time  has  been  varied  from  .05  T  to  .  15  1  with 

CO  00 

negligible  effect  on  the  drag  calculation  as  shown  in  Figure  6.  The  effect 
of  this  change  on  density  and  temperature  on  the  stagnation  line  is  shown 
in  Figures  III-l  and  III  -  2 .  There  is  no  apparent  effect  on  density  but 


-fc;  distance  along  stagnation  line 


Figure  III -1  Effect  of  Calculation  Parameters  on  Stagnation  Line  Density 
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X^DISTANCE  ALONG  STAGNATION  LINE 


Figure  III- 2  Effect  of  Calculation  Parameters 
on  Stagnation  Line  Temperature 

temperature,  which  is  more  sensitive,  shows  a  trifle  more  scatter  at 
Atm  =  .15  tm. 

Cell  size  has  not  been  varied  systematically  so  that  its  effect  could 
be  completely  isolated  but  the  results  of  many  different  calculations  with 
different  cell  sizes  suggest  that  as  long  as  the  cell  width  is  not  greater 
than  two  mean  free  paths  in  regions  of  strong  gradients  the  results  are 

not  strongly  influenced. 

The  calculation  is  performed  by  instantaneously  inserting  the  body 
and  allowing  the  field  to  evolve  to  an  approximate  steady  state  at  Time 
=  T  .  The  calculation  of  the  cumulative  average  to  determine  flow 
properties  is  begun  at  this  point.  The  initial  number  of  molecules  in  each 
cell  thus  affects  the  sample  size  involved  in  determining  this  approximate 
steady  state,  and  hence  the  accuracy  of  the  calculation  at  some  later 
time,  since  dispersion  in  the  initial  field  takes  some  time  to  be  averaged 
out.  The  density  and  temperature  profiles  on  the  stagnation  line  of  a 
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cylinder  at  time  =  ig^,  for  25  and  9  molecule  s/cell  are  shown  in  Figures 
II1-3  and  III -4.  Tue  larger  sample  has  reduced  the  initial  scatter  in 
density  but  temperature,  which  is  more  sensitive,  shows  large  scatter 
for  both  values. 

The  effect  of  field  size  is  shown  in  Figures  III-5  and  III-6  where 
density  and  temperature  profiles  on  the  stagnation  line  and  along  a  line 
which  coincides  with  the  upper  boundary  of  the  small  field  are  compared. 
The  field  sizes  are  shown  in  Figure  8.  There  is  negligible  effect  on  the 
stagnation  line  in  front  of  the  body  and  only  a  very  slight  effect  appearing 
on  the  downstream  portion  of  the  upper  boundary  line.  Temperature  does 
appear  to  be  more  sensitive  than  density  and  there  is  some  unexplained 
scatter  in  the  temperature  values  on  the  stagnation  line  in  the  rear  of  the 
body.  The  upper  boundary  condition  is  approximate  and  molecules  cross¬ 
ing  the  boundary  into  the  field  would  be  expected  to  introduce  something 
extraneous  to  the  disturbance  from  the  body;  however,  since  the  flow  is 
hypersonic  or  strongly  supersonic,  the  effects  of  errors  introduced  by  the 
feeble  thermal  motion  across  the  boundary  are  swept  downstream  rapidly 
and  show  up,  if  at  all,  only  near  the  boundary  downstream  of  the  body. 


FigureIII-3  Initial  Field  Scatter  Density  on  Stagnation  Line  of  a  Cylinder 
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Figure  III-4  Initial  Field  Scatter  Temperature  on  Stagnation 
Line  of  a  Cylinder 
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Figure  III-5  Effect  of  Field  Size  on  Density  in  the  Flow  Field  of  a  Cylinder 
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Figure  III- 6  Effect  of  Field  Size  oh  Temperature  in  the  Flow  Field 
of  a  Cylinder 
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APPENDIX  IV 


DESCRIPTION  OF  LOW  DENSITY  WIND  TUNNEL 
EXPERIMENTAL  SET-UP 

IV.  1  EXPERIMENTAL  METHODS  AND  PROCEDURES 
IV.  I.  1  Experimental  Equipment 

All  experiments  were  carried  out  in  the  TRW  Systems  Low  Density 
Wind  Tunnel.  This  facility  is  a  steady  flow  tunnel  exhausted  by  two  diffu¬ 
sion  pumps  backed  by  a  mechanical  vacuum  pumi;.  Two  nozzles  were 
used  in  these  experiments;  one  has  an  area  ratio  of  twenty-five  with  an 
exit  diameter  of  tv  o  inches,  and  is  conical  with  a  ten  degree  half  angle. 

The  other  has  an  area  r.  Uo  of  188  with  an  exit  diameter  of  1.  7  inches  and 
is  also  conical  with  \  ten  degree  half  angle.  Both  argon  and  nitrogen 
were  used  as  the  working  gas,  and  were  supplied  at  room  temperature 
from  high  pressure  botties.  The  experiments  were  carried  out  0.  4  inches 
downstream  of  the  nozzle  exit  plane  in  the  large  low  pressure  test  section. 
The  flow  conditions  used  at  this  section  in  the  flow  may  be  found  in 
Tables  IV -1  and  IV- 2. 

IV.  1.  2  Drag  Measurements 

The  drag  balance  used  in  thi6  experiment  was  an  adaptation  of  the 
Cahn  Electrobalance,  Model  RG,  which  is  a  remote  control  unit,  capable 
of  operation  in  vacuum  and  has  a  very  high  sensitivity.  The  balance  is  a 
null-type  single  component  type  beam  balance.  The  beam  is  mounted  on 
a  meter  movement  which  can  be  driven  to  a  null  position.  The  null  posi¬ 
tion  is  determined  by  a  photo  cell  which  senses  light  striking  it  through  a 
slit.  In  front  of  the  slit  the  balance  beam  moves  up  and  down,  thus  allow¬ 
ing  more  or  less  light  to  reach  the  cell.  If  an  amount  of  light  other  than 
that  of  the  null  position  reaches  the  ceil,  a  feed  back  circuit  drives  the 
meter  movement  to  return  the  beam  to  its  null  pof  ition.  "lhe  current 
necessary  to  accomplish  this  return  has  been  calibrated  and  can  be  read 
out  directly  as  a  force. 

In  order  to  use  the  balance  in  this  experiment,  certain  modifications 
had  tt  be  made.  A  picture  of  the  balance,  see  Figure  IV-1,  will  help  to 
explain  these.  In  order  to  support  the  sting  necessary  to  mount  the  model, 


a  clamp  had  to  be  made  which  holds  the  sting  arrangement  to  the  balance 
beam.  The  sting  then  goes  through  a  shield  and  finally  attaches  to  the 
model.  In  order  to  minimize  the  weight  effects  of  the  clamp  and  sting  a 
counter  weight  was  built  which  balances  the  moments  about  the  pivot.  It 
is  also  necessary  to  lower  the  center  of  gravity  of  the  system  below  the 
pivot  in  order  to  be  operating  in  a  stable  condition.  To  support  the  bal¬ 
ance  in  the  tunnel,  a  case  was  built  which  could  be  moved  in  two  directions 
on  the  support  pad  and  on  which  the  sting  shield  could  be  mounted.  This 
sting  shield  is  necessary  to  keep  down  the  drag  forces  on  the  sting  as 
much  as  possible. 

In  Figure  IV-2  can  be  seen  the  actual  support  of  the  models  on  the 
balance.  The  model  is  supported  on  the  upright  part  of  the  sting  by  a 
very  thin  horizontal  piece.  In  order  to  investigate  the  size  effect  of  this 
horizontal  piece  t\.o  different  diameter  stings  were  used.  One  was  a  1.  5 
inch  long,  0.  028  inch  diameter  stainless  steel  tube;  the  other  a  1.  5  inch 
long,  0.  010  inch  diameter  tungsten  rod.  On  the  0.  125  inch  diameter 
models,  this  gave  a  sting  to  base  area  ratio  of  5%  and  0.  6%  respectively. 


Figure  IV- 1.  Balance  and  Model  in  Test  Configuration 
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In  order  to  eliminate  the  drag  of  the  support  mechanism  from  the 
model  drag  a  dummy  model  and  sting  arrangement  was  made.  Again  in 
Figure  IV-2  we  see  a  schen  tic  of  this  set-up.  A  dummy  sting  of  the 
proper  size  is  inserted  in  the  beam  clamp  and  positioned  in  the  flow.  A 
model  was  then  attached  to  an  external  support  free  from  the  balance  and 
positioned  in  front  of  the  dummy  sting.  In  order  to  give  correct  readings, 
the  model  was  put  within  0.  001  inches  from  the  sting.  Since  the  balance 
is  of  the  null  type,  the  sting  always  remains  in  the  same  position  even 
though  varying  forces  are  applied  to  it  during  a  run.  The  effect  of  the 
sting  holding  the  dummy  model  was  determined  to  be  negligible 

Two  models  were  used  in  these  drag  investigations.  One  was  a 
0.  125  inch  diameter  nylon  sphere,  the  other  a  9°  half  angle  aluminum 

cone  with  a  base  diameter  of  0.  125  inches  and  a  nose  bluntness  ratio 
of  0.10. 

IV .1.3  Flow  Field  Surveys 

The  flow  field  in  the  stagnation  region  of  the  various  models  tested 
was  investigated  with  a  free  molecule  equilibrium  temperature  probe. 
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Schematic  of  Model  and  Dummy  Model  Setup  in  the  Balance 
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Figure  IV-2. 


(References  4  and  5.)  The  models  used  in  this  survey  were  two  wedges  and 
four  cylinders.  One  wedge  was  sharp  with  a  10°  half  angle;  the  other  a  10 
percent  nose  bluntness  wedge  with  9°  half  angle.  Both  of  these  wedges 
were  run  at  zero  degrees  angle  of  attack. 

The  cylinders  had  diameters  of  0.  002,  0.  005,  0.  0153,  and  0.  125 
inches. 

The  free  molecule  probe  was  made  in  the  form  shown  in  Figure  IV-3. 
A  small  current,  about  0.  08  ma  was  passed  through  the  vertical  tungsten 
wire  (0.  0002  inch  in  diameter).  The  voltage  drop  in  the  central  0.  5  inches 
of  the  wire  was  measured  through  potential  leads  made  of  0.  00025  incti 
diameter  platinum-tungsten  wire.  The  voltage  drop,  or  resistance,  which 
changes  with  temperature  could  then  be  recorded.  In  addition  to  providing 
the  potential  leads,  the  platinum- tungsten  wires  form  thermocouples  with 
the  tungsten  wire.  The  two  junctions  could,  with  the  current  shut  off  in 
the  central  wire,  be  used  to  record  the  point  temperature  at  the  ends  of 
the  resistance  section.  The  probe  was  calibrated  in  an  oven  before  using. 
The  center  part  of  the  wire,  including  the  two  thermocouples,  was  inserted 


Figure  IV-3.  Schematic  of  the  Equilibrium  Temperature 
Probe  and  Circuitry 


into  a  small  copper  block  which  was  heated  by  an  electric  resistance 
heater.  This  was  done  to  simulate  the  conditions  existing  while  testing  in 
the  wind  tunnel.  The  wire  output,  both  the  Emf  and  thermocouple  read¬ 
ings,  were  then  plotted  on  an  x-y  recorder  against  a  Cr-A^  thermocouple. 
A  calibration  curve  for  the  wire  Emf  and  thermocouple  output  against 
temperature  was  then  constructed. 

The  probe  was  made  by  stretching  all  three  wires  over  the  po6ts  and 
springs  (Figure  IV-4)  and  securing  them  under  tensions  of  approximately 
two-thirds  of  their  yield  stress.  Such  tensions  are  necessary  to  minimize 
deflections  under  drag  loads.  The  junctions  are  then  copper -plated  and 
soft  soldered  to  form  the  joints.  Care  was  taken  to  make  the  points  as 
small  as  possible  6ince  the  recovery  temperature  of  the  system  is  depend¬ 
ent  on  the  Knudsen  number.  Figure  IV-4  shows  the  probe  in  position 
between  the  model  and  nozzle  exit  plane. 

In  the  experiment  the  wire  probe  was  generally  held  stationary  and 
the  models  moved  on  a  lathe  cross  feed.  The  whole  stagnation  region 
could  thus  be  investigated.  The  model  location  was  monitored  by  d5al 


Figure  IV-4.  Equilibrium  lemperature  Probe  in  Operating 
Condition  in  Wind  Tunnel 
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indicators  of  0.  0C1  inch  least  count.  The  x-axis,  or  stagnation  line, 
location  of  the  model  was  also  indicated  by  a  Sanborn  displacement  trans- 
ducer  whose  output  was  directly  connected  to  the  x-axis  of  an  x  -  y 
recorder.  The  resistance  wire  and  thermocouple  outputs  were  connected 
to  the  y-axis  of  the  recorder. 

In  the  case  of  the  blunt  9°  wedge  a  modification  had  to  be  made. 

Since  this  model  was  run  both  at  liquid  nitrogen  temperature  in  addition  to 
the  recovery  temperature  run,  it  was  necessary  to  hold  the  model  as 
firmly  as  possible  to  avoid  deflection  under  liquid  nitrogen  conditions.  To 
accomplish  this  it  was  necessary  to  mount  the  wire  probe  on  the  crossfeed 
and  hold  the  model  firmly  to  the  tank.  The  wire  position  was  monitored  in 
the  same  way  the  model  had  been  in  the  other  runs. 

'’’he  temperature  in  the  reservoir  was  measured  by  a  chromel-alumel 
thermocouple,  and  the  model  and  wire  holder  post  (which  forms  the  cold 
junction  of  the  wire  thermocouples)  were  measured  by  copper -constantan 
thermocouples. 

The  0.  125  inch  cylinder  was  made  from  a  steel  rod  and  supported 
much  like  the  model  in  Figure  IV- 4.  The  smaller  cylinders  were  actually 
wires  supported  by  a  U-shaped  frame  attached  to  the  cross  feed. 

The  10°  sharp  wedge  was  made  from  oil  hardened  steel  and  ground 
to  a  sharp  edge  of  approximately  0.  0001  to  0.  0003  inches.  The  blur* 
wedge  was  also  machined  from  steel  with  a  hole  drilled  through  the  back. 
The  liquid  nitrogen,  supplied  from  a  standard  pressurized  bottle  and 
vented  to  atmosphere,  was  passed  through  that  hole. 

The  alignment  of  model  and  wire  was  accomplished  by  observing  the 
adjustments  through  a  jig  transit  with  an  optical  micrometer.  Wire  and 
model  were  aligned  parallel  within  0.  001  inch  over  the  length  of  the  model. 
The  jig  transit  was  also  used  to  measure  the  deflection  of  the  wire  under 
drag  loads.  Some  typical  deflection  values  are  presented  in  Table  IV-5. 

Probe  deflection  becomes  a  problem  when  one  moves  off  the  stag¬ 
nation  line  since  the  drag  on  the  probe  wires  changes  as  one  moves  to 
different  parts  of  the  flow  field.  Since  it  would  be  a  very  time-consuming 
process  to  measure  the  wire  deflection  at  all  the  conditions  used  in  the 
experiment,  it  was  decided  to  use  the  information  gathered  on  the 
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stagnation  line  surveys  to  reduce  off  line  data.  The  correction  was 
assumed  to  be  constant  for  equal  densities  or  temperature  ratios,  and  was 
carried  out  accordingly.  Spot  checks  of  this  method  showed  it  to  be  satis¬ 
factory  and  the  method  was  used  throughout  the  experiments. 

IV.  2  DATA  REDUCTION 

IV.  2.  1  Flow  Field  Calculations 

The  usual  method  of  flow  field  calculation,  using  the  impact  and 
stagnation  pressure  measurements  with  an  assumption  of  isentropic  flow, 
was  used  in  reducing  the  data.  Surveys  were  run  in  the  test  region  of  the 
nozzle,  both  axially  and  radially.  Since  the  nozzles  were  conical  and  not 
contoured,  both  radial  and  axial  Mach  gradients  exist.  A  range  of  stag¬ 
nation  pressures  was  found  for  both  nozzles  where  the  flow  field  consisted 
of  a  large  enough  isentropic  core  to  envelop  the  free  molecule  pro’  and 
where  the  gradients  were  small  over  distances  comparable  to  the  model 
dimensions.  The  calculations  were  then  carried  out  using  the  tables  of 
Reference  15. 

The  viscosity  of  air  was  determined  by  the  Sutherland  formula 

T  3/2 

=  2'  27  T.  ”"~9TT  x  10’8  lb  Sec/f*2 


For  Argon  the  viscosity  was  interpolated  from  the  tables  in  Reference  16. 
The  mean  free  path  was  then  found  from  the  equation 


\ 


03 
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The  flow  field  parameters  which  are  of  importance,  and  against 
which  the  data  were  plotted  are  as  follows. 


For  the  drag  measurements  the  drag  coefficients  were  plotted 
against  Reynolds  number.  In  the  case  of  the  sphere  drag  the  Reynolds 
number  was  based  on  free  stream  values  and  is  defined  as 


VpD 
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For  the  cone  drag  the  Reynolds  number  suggested  by  Potter  (Reference  17) 
was  used: 


or 


where 


1/2 


Re 


we 


w 


where  i  is  the  wetted  length  of  the  cone  and  is  evaluated  at  the  wall 
temperature.  This  Reynolds  number  was  chosen  because  it  seems  to 
correlate  all  cone  data  in  the  transition  regime,  independent  of  bluntness 
or  wall  to  stagnation  temperature  ratio. 

IV.  2.  2  Drag  Coefficient  Determination 

Since  the  setup  for  both  the  spheres  and  cones  was  the  same,  the 
data  were  reduced  identically.  Force  at  the  pivot  point  is  read  off  the 
control  box  of  the  micro  balance.  To  obtain  the  force  acting  on  the  model 
the  proper  moment  arms  have  to  be  used.  The  sting  correction  numbers 
are  obtained  in  a  similar  fashion.  Drag  coefficient  is  calculated  in  the 
following  way: 

F  -  F 

r  model  +  sting  sting 

D  "  "qA 

where  q  is  the  dynamic  pressure  and  A  the  area  based  on  the  model 
diameter. 
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Table  IV-3.  Sample  Deflection  Measure  aents  for  the 
Equilibrirm  Temperature  Probe 


Model 

Gas 

Stagnation 
Pressure 
Pt,  mmHg 

Mach 

Number 

M 

Distance 
From 
Stagnation 
Point  xf 
inches 

Deflection 
at  x 
inches 

0.  0156  inch 

N2 

100 

6.  35 

0.  040 

Cylinder 

0.  030 

0.  020 

0.  010 

0.  002 

0.  005 

0.  001 

250 

6.  75 

0.  040 

1 

0.  030 

1 

0.  020 

0.  010 

0.  005 

Ar 

15 

5.  20 

0.  040 

0.  004 

0.  030 

0.  020 

0.  010 

1 

0.  005 

25 

5.43 

0.  040 

0.  030 

0.  020 

0.  005 

0.  010 

0.  002 

0.  005 

0.  003 

40 

5.  72 

0.  040 

0.  030 

-  - 1 

0.  020 

1 

0.  010 

0.  005 
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